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The general setting for Shape Analysis 

Sylvain Arguillere* 


Abstract 

In shape analysis, the concept of shape spaces has always been vague, requiring a case-by¬ 
case approach for every new type of shape. In this paper, we give a general definition for an 
abstract space of shapes in a manifold M. This notion encompasses every shape space studied 
so far in the literature, and offers a rigorous framework for several possible generalizations. 
We then give the appropriate setting for LDDMM methods of shape analysis, which arises 
naturally as a sub-Riemannian structure on a shape space. This structure is deduced from 
the space of infinitesimal deformations and their infinitesimal action. We then describe the 
properties of the Hamiltonian geodesic flow, and study several applications of equivariant 
mappings between shape spaces. 


Introduction 

Mathematical shape analysis is a relatively recent area of study, whose goal is to compare several 
shapes while keeping track of their geometric properties. This is done by finding a deforming an 
initial shape (a template ) into another target shape in a way that minimizes a certain cost that 
depends on the properties of the shapes. This implies that a cost has been assigned to every 
possible deformation of the template, the design of this cost function being a crucial step in the 
method. This technique has lead to great contributions in computational anatomy (see [29]). 

A powerful method, the Large Deformation Diffeomorphic Metric Mapping (LDDMM), rep¬ 
resents deformations as flows of diffeomorphisms on the ambient space M in which the shape is 
embedded, this flow being generated by certain time-dependent vector fields [25, 49, 50] on the 
ambient space. These vector fields are chosen at each time in a fixed space V of ’’infinitesimal 
transformations”, equipped with a Hilbert product (•,•). In other words, we have a deformation 
t M- q(t) of the initial shape g(0) = qo given by 

q(t) = <p(t) ■ q 0 , 

where tp(t) is the flow at time t of the time-dependent vector field ( 4 1(1) € F: dt<p(t,x) = 
X(t,<p(t,x)). The cost of this deformation is then defined by the energy, or action, of the defor¬ 
mation given by 

J (X(t),X(t))dt. 

The goal is then to get the initial shape close to the target shape while minimizing this cost. Note 
that these transformations preserve local (such as the smoothness) and global (such as the number 
of self-intersections) geometric properties of the deformed shape. 

One of the main issues with mathematical shape analysis is that the concept of a “shape space” 
is rather vague: the literature only studies examples of shape spaces (usually spaces of embeddings 
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of a manifold N in M, or just spaces of submanifolds in M [13, 52]), without giving a good definition 
of what a shape space is. This not only forces a case by case approach, it can also limit the kind 
of shapes that are studied. For example, none of the existing methods tackle the issue of fibered 
shapes (submanifolds of M coupled with a fixed tangent vector field, representing the direction of 
a fiber), which would be very useful when studying the movement of a muscle, for example. This 
paper aims to remedy this situation by giving an abstract definition of a shape space that unifies 
and generalizes all the different shape spaces studied so far, and to generalize to this setting all 
the tools of the LDDMM framework. 

In order to define the space of infinitesimal transformations and the corresponding norm, two 
ways are classically considered in the existing literature. The first, used for example in [13, 11, 39], 
consists of considering all smooth vector fields with compact support. The resulting structure is 
a weak right-invariant Riemannian structure on the Lie group of smooth diffeomorphisms with 
compact support on M [35], which projects into a weak Riemannian metric on the shape space. 

The second, which is the point of view adopted throughout most of this paper, consists of 
considering an (a priori arbitrary) Hilbert space (V, (■,•)) of vector fields on M with continuous 
inclusion in the space of vector fields vector fiels of high enough Sobolev regularity. The space V is 
usually defined using a corresponding reproducing kernel (see Section 1.2). This is the framework 
that lead to the development of the so-called LDDMM methods (see [8, 9, 10, 14, 20, 33, 40, 41, 53]). 
The shape spaces are then usually studied as strong Riemannian manifolds , for which deformations 
with minimal costs are, again, geodesics. 

However, we will see that in this context, the shape spaces were actally sub-Riemannian, and 
not Riemannian. In the existing literature, the authors were not familiar with sub-Riemannian 
geometry, which is why the subject was treated from the Riemannian point of view 1 . Moreover, 
since the shape spaces we are dealing with are infinite dimensional, we get infinite-dimensional 
sub-Riemannian manifolds. This complicates matters even more, but can be worked around, as 
we will see. 

The paper is organized as follows. In Section 1, we recall the results of [5] on right-invariant 
sub-Riemannian structures on the group Diff(M) of diffeomorphisms of a smooth manifold M with 
bounded geometry, which will be necessary for all the results of this paper. Then, in Section 2, we 
give the definition of an abstract shape space in M, deduce the sub-Riemannian structure it inherits 
from the action of Diff(M), and study the corresponding notions of distance, geodesics, singular 
curves, and their applications to LDDMM methods. We also give some examples of Hamiltonian 
geodesic equations for several shape spaces, some that were already known but treated from a 
Riemannian viewpoint, some completely new. Finally, in Section 3, we look for symmetries in the 
geodesic flow, and study various applications of equivariant mappings between shape spaces. We 
also define the notion of lifted shape spaces which generalize the concept of diffeons from [54] , and 
give a few examples of geodesic equations on finite dimensional lifted shape spaces. 

1 Sub-Riemannian structures on groups of diffeomorphisms 

In shape analysis, and in the LDDMM methods in particular, one uses flows of vector fields 
on the ambient space M to compare shapes. Moreover, these vector fields belong to a fixed, 
arbitrary Hilbert space V. However, as shown in [5], such a setting naturally endowes the group of 
diffeomorphisms of M with a sub-Riemannian structure. The purpose of this section is to describe 
this structure and to give a brief summary of the results of [5]. 

1 This does not make their result incorrect though: first of all, in these methods, geodesics are usually found using 
the Hamiltonian viewpoint. Second of all, in practical applications, the shape spaces are often finite dimensional 
spaces of landmarks, for which the structure is actually Riemannian. 
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1.1 The group of diffeomorphisms 

We consider (M, g M ) and (IV, g N ) smooth oriented Riemannian manifolds with bounded geometry. 
Recall that this means that their global injectivity radius is positive and that for every ieN, the 
Riemannian norm of the *-th covariant derivative of their curvature tensor is bounded. Under 
these assumptions, for any integer s > dim(M)/2, one can define a smooth Hilbert structure on 
the space H S {M,N) of maps of Sobolev class H s from M to N (see [5, 26, 47]). These spaces are 
metrisable and we denote (m,n) the corresponding distances. 

Remark 1. We obtain a decreasing sequence of Hilbert manifolds with dense continuous inclusions 
H S+1 (M, N) H S (M,N), so that the limit H°°(M,N) is a smooth Inverse Limit Hilbert (ILH) 
manifold [43]. 

For s > d/2 + 1, let V S (M) = (H S (M, M) flDiff(M)) e be the connected component of e = Mm 
in the space of C 1 -diffeomorphisms of M that also belong to H S (M, M). Then V s (M) is an 
open subset of H S {M,M), and hence a Hilbert manifold. Its tangent space at e is given by 
T e T> s (M) = V s (TM) the space of all vector fields on M of Sobolev class H s . It is given everywhere 
else by T v V s (M ) = T S (TM) oip. It is also a group for the composition (<p,ip) >-» ipotp. This group 
law satisfies the following conditions: 

1. Continuity: (<p,ip) i —> <p o ip is continuous. 

2. Smoothness on the left: For every ip £ T> S (M), the mapping : ip M- (p o ip is smooth. 

3. Smoothness on the right: For every k £ N \ {0}, the mappings 

V s+k (M) x D S (M) —> D S (M) T s+k (TM) x V s (M) —> TV S [M) . 

(<* 2 , VO 1 —> +°ip {X,ip) i —> Xoip 


are of class C k . 

Remark 2. In particular, the intersection of the sequence (2? s (M )) s>d / 2 +1 is equipped 

with an ILH Lie group structure, on which composition of diffeomorphisms is smooth (and the 
above conditions are still valid for s = oo) [27, 43, 47]. 

Curves on V S (M). Let s be an integer greater than d/2+ 1. For every tp(-) £ H 1 ( 0,1 
define the logarithmic velocity of </?(•) by 


x (-) = <f(-) ° <*2(0 _1 G L 2 (0, 1; T S (TM)). 

Note that, by construction, we have <p(t) = X(t) o <p(t) for almost every t £ [0,1]. 

In other words, any curve </?(■) £ iL 1 (0,1; V S (M)) of diffeomorphisms is the flow of a time- 
dependent right-invariant vector field on T> S (M) with square-integrable LC-norm. Conversely, 
thanks to the measurable version of the Cauchy-Lipschitz theorem (see, e.g., [48]), any time- 
dependent vector field X(-) £ L 2 ( 0,1;T S (TM)) generates a unique flow ip x (-) £ H 1 ( 0,1;27 S (M)) 
such that <^' X (0) = e, and therefore a unique curve ip(-) = ip x (•) o<^q for any fixed intitial condition 
ipo £ V 8 {M). A group that satisfies such properties is called regular [35]. In other words, V S {M) 
satisfies the additional property: 

4. Regularity: For any tpo £ T> S (M) and X(-) £ L 2 (T S (TM), there is a unique curve tp{-) £ 
H l ( 0,1 ;V S (M)) such that <^(0) = ipo and ip(t) = X(t) o ip(t) almost everywhere on [0,1]. 
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1.2 Sub-Riemannian structures on V s (M). 

Let us first define general sub-Riemannian structures on a Banach manifold. 

Definition 1 ([4]). Let At be a smooth Hilbert manifold. A strong sub-Riemannian structure of 
class C k on At, k £ N U {oo}, is a triple (£,g,£), where tt : £ —>- M is a smooth Hilbert vector 
bundle on M endowed with a smooth, strong (i.e. fiberwise Hilbert) Riemannian metric g, and 
£ : £ —> TM a vector bundle morphism of class C k . 

A horizontal system is a curve t £ [0, 1] (q(t), u(t)) £ £ of class L 2 , such that its projection 

t i —> q(t) to At is of Sobolev class H 1 and satisfies for almost every t £ [0,1], q(t) = £, q (t)u(t). Its 
length and action are respectively defined by 


L{q,u) = L A (<L U ) = \ J g q (t){u{t),u{t))dt. 


A horizontal curve is the projection q(-) to At of a horizontal system. 


Note that, if £ is one-to-one, there is a bijective correspondance between horizontal systems 
and horizontal curves. 


Definition 2. For At = T> S (M), we can define strong right-invariant structures as follows: fix 
V an arbitrary Hilbert space of vector fields with Hilbert product (•, •} and continuous inclusion in 
T s+k (TM), k £ N. The sub-Riemannian structure induced by V on V S (M) is 

SRy(s) = {V s (Af) x V, (-, -),dR), 


where, for (<p,X) £ V s (M) x V, 

dR v {e).X = dR v X = X o<p £T V D S {M). 

SRy(s) is a strong sub-Riemannian structure of class C k . See [5] for various examples. 

Remark 3. Let V* denote the dual space of V, and K v : V* — > V the inverse of the canonical 
isometry X £ V i —> (X , •) £ V*. If V has continuous inclusion in the space of continuous vector 
fields on M, then Ky is given by convolution with a reproducing kernel K. Such a kernel is a 
mapping (x, y) £ M 2 i —> K(x , y), with K(x, y) an linear morphism from T*M to T X M. It is defined 
by the property that for every (y,p) £ T*M, the vector field x ^ K(x,y)p belongs to V and 

VX£V, (K(;x)p,X)=p(X( X)). 

Then, if an element P of the dual V* can be represented by a one-form with (distributional) 
coefficients 2 so that we can write 


P(X) = I P{x){X{x))dx, 
Jm 


then the vector field Y £ V given by 


Y(x)= [ K(x, y)P(y)dy 
Jm 


satisfies 

VXeV, (Y,X) = P(X). 

2 For example, T S (TM)* = T~ a (T*M ) can be identified with the space of one-forms on M with coefficients in 

H~ S (M). But since we assume that V has continuous inclusion in PfTM), any such 1-form also belongs to V*. 
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Conversely, such a mapping K is the reproducing kernel of a unique Hilbert space of vector fields 
with continuous inclusion in the space of continuous vector fields if and only if it satisfies 


du(x)(K(x, y)dv{y)) = 0 =>■ Vx £ M, 


K(x,y)dv(y) = 0 


JJMxM J M 

for every compactly supported 1-form v with coefficients in the space of Radon measures on M 

[53]. 


Horizontal curves and end-point mapping. The flow ip x (-) of any X(-) £ L 2 ( 0, 1; V) is a 
horizontal curve for the structure defined in the previous paragraph, and conversely any horizontal 
curve ip(-) £ H 1 (0,1; V s (M)) with i^(0) = e is the flow of A(-) = ip(-)o £ L 2 (0,1; V). Hence, 
we can identify the set of all such horizontal curves with the Hilbert space L 2 (0,1; V) through the 
correspondance X o ip x . Note that the length and action of such a curve are given respectively 

by 

L(<p x (-))=L(X(-)) = f y/(X(t),X(t))dt, 

JO 

A(<p x (.)) = A(X(.)) = \j\x{t),X(t))dt. 

The endpoint mapping end : L 2 { 0,1; V) —> V S (M) is given by end(A(-)) = p x (1) and is of 
class C k (recall that k is such that V T s+k (TM) is a continuous inclusion). 


Sub-Riemannian distance. The sub-Riemannian distance dsR between two elements ipo and 
ip i of T> S (M) is defined as the infimum of the length of horizontal curves steering ipo to (pi, with 
the agreement that dsnipo^i) = +oo whenever there is no horizontal curve steering ipo to ip\. 
A curve y>(-) : [0,1] —>• V S (M) is said to be minimizing if dsii{<p(0), y>(l)) = L(*p{-)). As usual, a 
curve minimizing the energy among all other curves with the same endpoints also minimizes the 
length. 

Proposition 1. The sub-Riemannian distance d$R is indeed a distance (taking its values in 
[0, +oo]j, that is, dsRitpoiTi) = 0 implies <po = p\. Moreover, the corresponding topology on 
V s (M) if at least as fine as the intrinsic manifold topology. 

We also have the following generalization of Trouve’s result [49] . 

Theorem 1 (Completeness of the distance[5]). Any two elements ipo and ip\ of V s (M) with 
dsR^POiTi) < +°° can be connected by a minimizing horizontal curve, and ( V s (M),dsR) is a 
complete metric space. 


1.3 Geodesic flow on V S (M ) 

We keep the framework and notations used in the previous sections, but we assume k ^ 1. A 
geodesic < p x ( m ) starting from e is a critical point of the action A(A(-)) among all curves gfi (•) with 
the same endpoints. In other words, for every C 1 variation a G (—£,£) X a (-) G L 2 ( 0,1; V) such 

that end(A a (-)) = ip\, and Xq = X, we have 

d a (A(X a (-)) | Q=0 = 0. 
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The three kinds of geodesics. It easy to see that for any such curve, the couple of linear 
operators 

(dA(X(-)),dend(X(-))) : L 2 (0,l;f)^lx T Vl V s (M) 
is not onto. This can imply one of two possibilities: 

1. There exists (A, Pi) elx T*V S (M) \ {(0, 0)} such that 

dend(A(-))*Pi + AdA(X(-)) = 0. (2) 

This is a Lagrange multipliers relation, and can be used to give a Hamiltonian caracterization 
of ip x (•). We have two subcases: 

(a) Normal case: A ^ 0. In that case, we can take A = — 1, so that dA(A”(-)) = dend(W(-))*Pi. 
We say that (p x ( •) is a normal geodesic. 

Note that, conversely, for any X(-) £ L 2 ( 0,1; V), if we have dA(X(-)) = dend(X(-))*Pi 
for some Pi, then tp x (•) is indeed clearly a geodesic. 

(b) Abnormal case: A = 0. In that case, (2) implies that ip x (-) is a singular curve. We say 
that (p x (■) is an abnormal geodesic. 

2. The image of (dA(X(-)), dend(A(-))) is a proper dense subset of T tf>1 V s (M). There are no 
Lagrange multipliers. We say that <fi x (-) is an elusive geodesic. No Hamiltonian caracteriza¬ 
tion can be given. 

Note that this last case only happens because we work in infinite dimensions, and is a well- 
known problem in the community of control theory [31, 36, 37]. It is caused by a topological 
incompatibility between the manifold topology of V s (M) and that of the sub-Riemannian distance 
dsR- See [4, 5] for a more comprehensive discussion on the appearance of elusive geodesics and 
their consequences. Also see [22, 23, 15, 38, 46] for more on abnormal geodesics. 

Remark 4. When studying shape spaces, V is usually a dense subset of V s (TM). This implies 
that dend(A(-)) has at least dense for all X(-). In that case, there are no abnormal geodesics, 
only normal and elusive ones. 

Hamiltonian formulation. For A £ {0,1}, define the Hamiltonian of the motion P, X) = 

P(X o ip) — ~(X,X), with (tp,P) £ T*V S {M) and X £ V. Then one can prove the following [5] 
partial Pontryagin maximum principle [44] . 

Proposition 2. Let ip(-) be the flow of X(-) £ L 2 (0,1; V), with tp( 1) = ipi, and (A, Pi) £ {0,1} x 
T*i T> S (M) \ {(0,0)}. Let (</?(•), P(-)) be the unique lift of ip(-) to T*V S (M) that solves the linear 
Cauchy problem P( 1) = Pi and P(t) = —d v Hx(ip(t), P(t), X(t)) for almost every t £ [0,1]. Then 

AcL4(<p(-)) = dend(AT(-))*Pi •<=>> K v dR* v(t) P{t) = A X(t) a.e. t £ [0,1], (3) 

In this case, we call P(-) a singular covector when A = 0 and a normal covector when A = 1. 

In particular, a curve <p(-), flow of X(-), is singular if and only if it can be lifted to a curve 
(<£>(•), P(-)) such that, for almost every t £ [0,1], 

P(t) = —( dX(t ) o ip(t))*P(t), V C ker P[t). 

On the other hand, define the normal Hamiltonian H : T*T> S (M) — > R by 


H(<p,P) = H 1 (cp,P,K v dR* v P). 
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H is of class at least C k . We see that normal geodesics are the projections to V s (M) of solutions 
on T*T> S (M) of the so-called Hamiltonian geodesic equation (<p,P) = (dpH((p,P),—d v H(ip,P)). 
Note that, using the reproducing kernel I\ of V, we have 

H {V,P) = \ [ P(x)(K(q>(x),<p(y))P{y))dxdy, 

z JMxM 

so that we get the usual formulas 
d P H(ip,P)(x) = [ K(ip(x),ip(y))P(y)dxdy, 

JM 

d v H(<p,P)(x) = \ [ (d 1 K(ip{x),ip(y))P(y))* P(x)dy+ \ [ (d 2 K{tp{y), tp(x))P(x))* P{y)dy. 
z Jm 1 Jm 

Moreover, one easily checks that (<p, P) —>• ( dpH(ip , P), —d v H(ip, P)) is a vector field of class C fe_1 
on T*T> S (M). We obtain the following theorem. 

Theorem 2. If k 2, then the Hamiltonian geodesic flow is well-defined, of class C k ~ l , and 
global. In other words, for every ( tpo,Po ) £ T*T> S (M), there is a unique solution (</?(•), P(-)) : K —> 
T*T> S (M) to the Cauchy problem 

0(0), P(0)) = Oo,P 0 ), 0(t),P(t)) = {d P H{y{t),P{t)),-d v H(ip{t),P{t))) a.e. t£ [0,1]. 

This solution is of class C k in time and C k ~ 1 in the initial conditions (<po; Po)- Moreover, any subarc 
of this solution projects, up to re-parametrization, to a normal geodesic on V s (M). Conversely, 
any normal geodesic is the projection of such a solution. 


Momentum formulation. We define the momentum map p : T*V S {M) —> T S (TM)* = T S (T*M) 
by p(tp, P) = (d R V )*.P. This proposition is proven in [5], and is connected to the EPDiff equation 
[7, 26, 32, 39], 

Proposition 3. Assume that k ^ 1. Then ip(-) € iP( 0,1;2? S (M)) be either a normal geodesic or a 
singular curve, with logarithmic velocity X(•) G P 2 (0,1; V), and let P(-) be a corresponding normal 
or singular covector. We denote by p(t) = p(<p(t), P(t)) the corresponding momentum along the 
trajectory. Then, for almost every t £ [0,1], we have 


fl(t) = ad * X ( t )P(t) = -C X (t)p(t). 


Here, adx : r !+ 1 (T¥) — > T S (TM), with ad X Y = and C x the Lie derivative with respect 

to X. As a consequence, 

p(t) = 0), 


for every t £ [0,1]. 


In this paper, we are also interested in the converse. 

Proposition 4. We assume that k 1. Fix p(-) £ H 1 {0,1-,'D S (AI)) a horizontal curve with 
logarithmic velocity X(-) £ L 2 (0,1;V). Let A £ {0,1} and t £ [0,1] p(t) £ T S (TM)* be 
continuous such that for almost every t £ [0,1], Kyp(t) = A X(t) and p(t) = ad * X ( t )p(t). Then ip(-) 
is a normal geodesic if A = 1 and a singular curve if A = 0, with covector P(t) = (dP*^)) _1 tO)■ 
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Proof. We already have KydR^^Pft) = A X(t) for almost every t. We just need to prove that 
P(t) = - d ip H x (ip{t),P{t),X(t )) = ~{dX(t) o ip(t))*P{t). Let Y G r s+1 (TM). We have 

P(t){Y o <p(t)) = -^( P(t)(Y o <£>(*))) - P(t)(dY o tp{t).X(t) o <p(t)) 
at s- v - s 

= Y]) - fi(t)(dY.X(t)) = -p(t)(dX(t).Y) 

= — P(t)(dX(t ) o ip(t).Y o ip(t)) 

= —(dX(t) o </?(f))*P(f)(H o <£(*)). 

Consequently, P(t) = —(dX(t)oip(t))*P(t) onT s+1 (TM)oip(t), and therefore on all of T lp 'D s (M) = 
r s (TM) o ip(t) by density. □ 

2 Shape spaces 

In this section, we give the definition for abstract shape spaces in M which allow the application 
of the usual LDDMM methods. We want the basic examples of shape spaces, such as spaces of 
landmarks, or spaces of embeddings of a fixed submanifold N, to be included in this definition. 
A common factor between those examples is that they are acted on by V s (M), and this action 
satisfies properties that are very similar to the particular topological group structure of V s (M). 
This action, along with the choice of a Hilbert space V of vector fields of M, endowes the examples 
with a length structure (considered Riemannian structure in the literature, but which is often only 
a sub-Riemannian structure), and the problem of comparing shapes can then be seen as a search 
for geodesics with respect to this structure. This leads to a Hamiltonian equation for the geodesic 
flow, which in turn lead to various optimization methods [52, 53]. 

Therefore, we define shape spaces as Banach manifolds on which the group of diffeomorphisms 
of M acts in a way that is “compatible” with this particular group structure. This then lets us 
project the sub-Riemannian structure induced on T> S (M) by an arbitrary Hilbert space of vector 
fields V to a sub-Riemannian structure on the shape space itself, for which horizontal curves are 
exactly curves induced by flows of elements of V. We study the corresponding sub-Riemannian 
distance, and prove in particular the geodesic completeness of the structure. Then, we find the 
Hamiltonian geodesic flow and the singular curves on the shape space for this structure, and their 
connection with the Hamiltonian geodesic flow and the singular curves of V s (M) itself. This allows 
the generalization of the usual LDDMM methods and algorithms to this broader category of shape 
spaces. 

2.1 Definition 

Throughout the section, M is a smooth Riemannian manifold of dimension d and bounded geom¬ 
etry. Let so be the smallest integer such that So > d/2. A shape space in M is a Banach manifold 
acted upon by V s (M) for some s in a way that is compatible with its particular topological group 
structure. 

Definition 3. Let S be a Banach manifold and <gN\ {0}, and s = Sq + £■ Assume that V S (M) 
acts on S , according to the action 


V S (M) x S —>■ 5 

ip-q 


( 4 ) 


We say that S is a shape space of order £ in M if the following conditions are satisfied: 



1. Continuity: (ip, q) >—>■ p ■ q is continuous. 

2. Smoothness on the left: For every q £ S, the mapping R q : p i—>■ p ■ q is smooth. Its 
differential at e = Id R d is denoted t; q , and is called the infinitesimal action of V s (TM). 

3. Smoothness on the right: For every k £ N, the mappings 

V s+k {M)xS —» S £ : T s+k (TM) x S 

(tp,q) 1 —> p ■ q ( X,q ) 

are of class C k . 

4- Regularity: For every X(-) £ L 2 (0,1; T S (TM)) and qo £ S, there exists a unique curve 
q(-) = q x (-) € i? 1 (0,1:5) such that q x (0) = q o and q x (t) = £ q x( t )X(t) for almost every t in 
[ 0 , 1 ]- 

A an element q of S is called a state of the shape. 

Remark 5. One can define shape spaces of order I — 0 (see [6]), but the action must still be 
performed by T> So+1 (M). 

Remark 6. T~L s ° +e (M, M) is a shape space of order I for the action by composition on the left. 
So is V So+e (M) itself. 

Remark 7. While we do not consider fractional order for shape spaces or diffeomorphisms of 
fractional Sobolev regularity out of a desire to keep things more simple, all definitions and results 
can be straightforwardly extended to those cases. 

Remark 8. The following properties immediately follow from the definition. 

• For every p G V s+k (M) and X G T s+k (TM), the mappings q H> p o q and q i-> f q (X) are of 
class C k . Moreover, q i-> £ q is of class C k ~ 1 (see [43] for example). 

• A shape space of order £ is also a shape space of every order £' > £. 

Definition 4. We say that q G S has compact support if there exists a compact subset U of M 
such that R q : p >-» p ■ q is continuous with respect to the semi-distance d H s 0 +e^ UM ^ on V s (M). 

In other words, q has a compact support if p ■ q depends only on the restriction of p to a 
compact subset U of M. 

Example 1. Here are some examples of the most widely used shape spaces: 

1. Let N be a smooth compact Riemannian manifold, and cto be the smallest integer greater than 
dim(iV)/2. Then S = Emb a ° +e (N 1 M) and Embf(iV, M), the manifolds of all embeddings 
q : N —» M respectively of class H a ° +i and C^, £ G N, are shape spaces of order max(l,£). 
In this case, 2? s o+ max ( 1 ^)(M) acts on S by left composition p ■ q = p o q, and this action 
satisfies all the required properties of Definition 3 (see [6] for the proof), with infinitesimal 
action f q X = X o q. Every q G S has compact support. 

2. A particularly interesting case is obtained when dim(S') = 0. Then N = {a\,... ,a n } is 
simply a finite set. In that case, for any £, the shape space S = C e (S,M ) is identified with 
the space of n landmarks in M: 

Lmk n (M) = {(zi,..., x n ) G M n \ Xi ± Xj if i ± j}. 

For q = (xi,... ,x n ), the action of V S ° +1 (M) is given by p ■ q = (p(xi),... ,p(x n )). For a 
vector field X of class H s ° +1 on M, the infinitesimal action of X at q is given by £, q {X) = 
(X(xi),... ,X(x n )). Spaces of landmarks are actually spaces of order 0 (also see [6]). 
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3. Let 5i and S 2 be shape spaces of respective orders £\ and £ 2 • Then the space S\ x S 2 is a 
shape space of order max(£i, £ 2 ) for the diagonal action ip o (g 1; q 2 ) = (ip ■ qi, <p ■ g 2 ). 

4. Let S be a shape space of order £ £ N. Then TS is a shape space of order £ + 1 for the action 
of D So+e+1 (M) on T<Si defined by <p ■ (q,v) = (<p ■ q,d q (<p ■ q){y)). 

Remark 9. The first example above implies that Emb°° (S,M) = n^Emb“ 0+£ (<S, M) is an ILH- 
manifold, on which U°°(M) acts smoothly as an ILH Lie group [43]. Such a space can be called a 
shape space of infinite order. They will be briefly studied in Section 3.2. 

Deformation of a state. A deformation of a state go is a curve g(-) : [0,1] —> S of the form 

Vie [0,1], q(t) = <p(t) ■ q 0 , 

where ip(-) : [0,1] — > V S (M) and g?(0) = e. Such a ip(-) is called a lift of the deformation q(-). Now 
if tp(-) is the flow of some vector field X(-) £ t 2 (0,l;T s (TM)), the deformation q(-) = ip(-) o q 0 
belongs to fL 1 (0,1;<S) and is the unique solution to the Cauchy problem 

Q{i)=€q(t)X(t)=X(t)oq(t), q(0) = q o , (6) 

for almost every t £ [0,1]. In other words, using the notations of 4 in Definition 3, for every 
X <EL 2 (0,1;V), and t £ [0,1], 

q x (t) = <p x (t) ■ q 0 . 

2.2 Sub-Riemannian structure on shape spaces 

Let S be a shape space of order £ ^ 1 in M, and fix k £ N \ {0}. Consider (V, (-, ■)) an arbitrary 
Hilbert space of vector fields with continuous inclusion in T So+e+k (TM). According to Section 1.2, 
we obtain a strong right-invariant sub-Riemannian structure SRy(s) of class C k induced by V on 
V s (M), with s = so + £■ 

The framework of shape and image matching. The classical LDDMM algorithms for exact 
shape matching seek to minimize 

(X(t),X(t))dt 

over every X £ L 2 (0,1;E) such that ip x ( 1) • qo = qi, where the template qo and the target q\ 
are fixed. Usually, one only wants to get ’’close” to the target shape, which is accomplished by 
minimizing 

(X(t),X(t))dt + g(ip x ( 1) • go) 

over every X £ L 2 (0, 1; V), where the endpoint constraint has been replaced with the addition of 
a data attachment term g(<p x (l) ■ qo) in the functionnal [28, 21]. The function g is usually such 
that it reaches its minimum at q±. 

The sub-Riemannian structure. The previous discussion leads us to define a sub-Riemannian 
structure on S as follows. 

Definition 5. The strong sub-Riemannian structure induced by V on S is the triple 

SR$ = (Sx V, <■,.), 0> 
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where f, in a slight abuse of notation, denotes the restriction to S x V of the infinitesimal action 
f : S x T S (TM) —>• TS. This structure is of class C k . Horizontal systems are couples 

{q{-),X{-))=H 1 {0,l-,S)xL 2 {0,l-,V) 

such that for almost every t £ [0,11, 

m= 

The curve q(-) is called a horizontal deformation of q{ 0). 

Remark 10. If f q (V) = T q S for every q £ S, this is actually a Riemannian structure. This is often 
the case in numerical simulations, where S is finite dimensional (usually a space of landmarks). 
However, in the general case, we do not obtain a Riemannian structure. 

The length and action of a horizontal system (g(-),X(-)) only depend on X(-) and are respec¬ 
tively given by 

L(q(-),X(-))=L(X(-))= [ y/(X(t),X(t))dt, 

Jo 

A(q(.),X(.)) = A(X(.)) = l -j\x(t),X(t))dt. 

The LDDMM algorithm can therefore be formulated as a search for sub-Riemannian geodesics on 
S for the structure SRy. 

A horizontal system (<;(•), Af(-)) satisfies q(-) = ip x (•) • q( 0). Moreover, is horizontal for 

SRy(s) on V 8 (M). The flow ip x (•) is called a horizontal lift of the curve q(-). This explains the 
term horizontal deformation for q(-). 

2.3 Minimal lift of a horizontal curve. 

The length and action of a horizontal system (q(-),X(-)) £ fly coincide with the sub-Riemannian 
length and action of the horizontal lift (p x (•) on V S (M). However, there may be two distinct 
controls Xi(-) ^ A^-) inducing the same horizontal curve q(-). For example, if S = Lmk n (M), 
and q 0 = {x \,..., x n ) £ S , then the flow of any X(-) £ L 2 (0,1; V) \ {0} such that 

X{t,Xi)= 0, i = 1,... ,n, t £ [0,1], 

leaves qo fixed, and induces the constant deformation t i —> q^. Since different controls can have 
different actions, it is natural to ask whether, given a horizontal deformation q(-), there exists a 
control Al(-) with minimal action such that q(-) = ip x (•) • g(0). 

To answer this question, fix q £ S. Then has a closed null set ker£ 9 C V. V is Hilbert, so 
ker f q admits a unique orthogonal supplement (ker^) 2 - such that the restriction of f q to (ker^j) 2 - 
is a bijection onto f q (V). The inverse £“ 2 : f q (V) (ker^) 2 * of this bijection, denoted in 
a slight abuse of notation, is the pseudoinverse of £ g . This means that = id^ 5 (v), but that 

is the orthogonal projection V —> (ker^) 2 -. 

Now, for every w £ f q V and X £ V with f; q V = w, we have 

(X,X)= min «y,y)) ^ x=ff l w. (7) 

YeV, ( g Y=w 

We immediately obtain the following lemma. 

Lemma 1. A horizontal curve q(-) £ H 1 (0,1; S) admits a unique control X(-) with minimal action 
A(X(-)) such that q(t ) = ^ q ^X(t), a.e. t £ [0,1]. This control is given for almost every t £ [0,1] 
by 
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The corresponding horizontal lift (p x (•) is called the minimal lift of q(-), and it is the horizontal 
curve on T> S (M) with minimal action such that 

q(-) = p X (-)-q(0). 

Remark 11. One can then define the length and action of a horizontal curve q(-) independently 
of a corresponding control by 

J 0 

However, L s and A 5 are much harder to handle than A and L (for example, they may not be 
differentiable), so we will not use them. 

The dual viewpoint. Computing the minimal lift of a given horizontal curve q(-) is very hard 
to do, as we would need to compute the pseudo-inverse £“ 1 . However, one can use a dual point of 
view to directly generate controls X(-) that correspond to minimal lifts. 

Lemma 2. Let q G M and p G T*M. The momentum map of the action ofV s (TM) on S is 
defined by p s : ( q,p ) H- £*p = poq G T S (TM)*. 

Then X = Ky£*p belongs to (ker^)- 1 . In particular, 

X = argmin{(y, Y) | Y G V, £ q Y = f q X). 

Proof. Let X = Kvfq'P (recall that Ky is the canonical isometry V* —> V). For every X 0 G ker(^ g ), 
we have (X, X 0 ) = (K^*p, X 0 ) = f,*p{X 0 ) = p(£ q X 0 ) = 0, hence X G ker^,)- 1 . □ 

We define the operator 


K q = £ q Ky? q : T^s c T q S. 

K q is symmetric (i.e., p\(K q p 2 ) = p 2 (K q pi)), positive semi-definite (i.e., p(K q p) ^ 0), and we have 
(Ky£*p, Ky£*p) = p(K q p). The map (q,p) i—>■ K q p is of class C k . Moreover ^ 9 (V) = K q {T*S) 
when £ q (V) is closed in T q S [18]. 

Then curves q(-) G H 1 { 0,1;5) satisfying 

q(t) = K q{t) p(t), p(t) G T* {t) S, 

almost everywhere are horizontal, and have minimal lift X(t) = Ky£*, t -.p(t), with action 

(X{t),X(t))dt = p{t) [K q(t) p{t)) dt. 

Examples using the reproducing kernel. As mentioned in Remark 3, if V is a Hilbert space 
of vector fields on M with continuous inclusion in V s {TM), then V is uniquely characterized by 
its reproducing kernel K, which is a section of the vector bundle L (T*M,TM) —» M x M such 
that J MxM P{x) (K(x,y)P(y)) dydx ^ 0 for every continuous one-form P : M —» T*M. More 
generally, this kernel also converts any continuous linear form P on V into the vector field X € V 
such that P = (X, •), by convolution of P with K: X(x) = K{x , -)P = f M K(x, y)P(y) d y. 

In shape analysis, it is common to start with an explicit kernel K instead of the space V. 
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Example 2. A commonly used example is the Gaussian kernel when M = W 1 : 

| u; _ y | 2 

K(x,y)p = e ^~p T , 

where we identified linear forms p £ (R d )* with lign matrices, and vectors in R rf with column 
matrices. 


Such a kernel does generate a unique Hilbert space of vector fields V, but the space itself is not 
explicit. Consequently, it is hard to say if a given vector field belongs to V. On the other hand, 
it is easy to generate explicit examples of elements of V : since we have an explicit formula for 
the reproducing kernel, we can simply take convolutions of continuous 1-forms (or, more generally, 
1-cocurrents) with K. This is another reason why the dual viewpoint is so important. 

Example 3. Take S = Lmk"(M). Any element q of S is a n-uple q = (aq,..., x n ) with aq £ M, 
and the infinitesimal action of I £ T S (TM) is A o q = £ q X = (X(aq),..., X(x n )). Any one- 
form on S is a n-tuple p = (ai,... ,a n ) with oq £ T*.M, and p(w) = Y2i=i a i( w i)- We have 
p o q(X) = p(X o q) = a i{X( x i)), and hence p o q = i a i ® d Xi , with S x the usual Dirac 
mass. Therefore, if (x,y) >->■ K(x,y) is the reproducing kernel of the Hilbert space V, we get 
K v£qP(x) = E"=i K{x, x i)ai, and 

n 

KqP = iqKvCqP = ^ {K{ X U X i) a i, * * * G T X± M X • • • X T Xn M. 

i =1 

In particular, 

n 

{ K vC q P, KvC q v) =p{K q p) = Pi( K ( X i,Xj)Pj). 

i,j =1 

Example 4. Take S = C°(N,M), with N a compact Riemannian manifold. An element q of S is 
a continuous mapping q : N —>■ M, and 

T q S = T 0 (q*TM) = K C°(N, TM) | Va e S, v{a ) £ T q(a) M} 

The infinitesimal action of X £ V at q is 


i q X = X oq£To(q*TM). 

A one-form p £ T*S is a section of q*T*M with (distributional) coefficients in the space of Radon 
measures. Then 


Now for X £ V, we have 


so for x £ M, 


For K q p : N —» TM, we obtain 


p(v) = / dp{a){v(a)). 

Jn 

P(tq X ) = f dp(a)(X(q(a))), 
Jn 

K v£, q p(x) = [ K(x,q(a))dp(a). 
Jn 


K q p(a)= / K(q(a),q(a'))dp(a') 


in 


and 


{KvCqP, KvCqP) = P{K q p) = [[ dp (a) (K (q(a ), q(a')))dp(a'). 

J Jn 2 
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2.4 Sub-Riemannian distance 


We keep the same notations as in the previous section. The sub-Riemannian distance induced 
by the sub-Riemannian structure SRy = (S x V, (•, ■),£) is denoted d RR : S x S —>• R+ U {oo}. 
Recall that dg R (qo,qi) is the infimum over the lengths of every horizontal system (q(-), X(-)) with 
g( 0 ) = 9 o and 9 ( 1 ) = 9 i- It is clear that d^ R is at least a pseudo-distance (it satisfies all axioms of 
a distance function except, possibly, the separation axiom). 

Note that we have 

d ffl( 9 o, 9 i)= inf (d SR (e,<p)) = d S R{e,R~^({qi})), ( 8 ) 

< p, <p-qo=qi 

with inf(0) = +00 and Rf^dqi}) = {<p £ V S {M) \ ip ■ q 0 = < 7 !}. Indeed, if dsR(e,ip ) < +00 
and (p ■ go = qi, we can consider X(-) such that p X (■) is the minimizing geodesic between e 
and p. But then dsR.(e,<p) = L(X(-)) ^ d RR (qo,qi). On the other hand, if (X n (-)) n6 N is a 
sequence of controls steering qo to q\ with L(X n (-)) —>• dg R (qo,q\) as n goes to infinity, then 
d-SR^Rgo 1 ^})) < d S R(e,p Xn ( 1)) ^ L(X n (-)). 

A horizontal system (q(-),X(-)) is said to be minimizing if it minimizes the action among all 
other systems with the same endpoints. The horizontal curve q(-) is called a minimal curve , while 
X(-) is a minimal control. In this case, Cauchy-Schwarz’s theorem tells us that 

df fl («(0),g(l)) = L(q(.),X(.)) = y/2A(q{.),X(.)). 

Note that the minimal lift </?(•) of a minimizing curve g(-) is a minimizing curve between e and 
<p(l). In this case, d§ R (q{0),q(l)) = d SR (e,p( 1)). 

Theorem 3. The sub-Riemannian distance is a true distance, that is, d RR (qo,qi ) = 0 implies 
do — qi- R induces a topology that is at least as fine as the manifold topology of S. 

Moreover, if qo £ S has compact support, and if we denote S qo = V s (M) ■ qo, then {S qo ,dg R ) 
is geodesically complete (any two points within finite distance can be connected by a minimizing 
geodesic). 

Proof. Since Banach manifolds are Hausdorff spaces, proving that d RR induces a finer topology on 
S will also prove that it is a true distance. For this, let U C S be an open subset for the manifold 
topology. Fix q £ 17. Then 7 ?” 1 (t/) is open in V s (M) and contains e. Theorem 1 implies that there 
exists r > 0 such that the sub-Riemannian ball B in V S (M) of center e included in I?“ 1 (C7). But 
we saw earlier in the section that R q (B) is exactly equal to the sub-Riemannian ball on S of center 
q with same radius as B. Therefore, U contains an open ball for d RRl and the sub-Riemannian 
topology is finer than the manifold topology. 

Now let us prove the geodesic completeness of our space for compactly supported states. Assume 
qo has compact support. It is clear that any element of S qo also has compact support. Hence, to 
prove that S qo is geodesically complete, we just need to prove that if d s (qo, qi) < + 00 , then qo and 
qi can be joined by a minimizing geodesic. 

Let (<7 n (-))neN be a minimizing sequence of horizontal curves in H 1 ( 0,1; S) such that q n ( 0) = qo 
and q n ( 1) = qi for every n, induced by a sequence of controls X n (-), such that \J2A(q n (-), X n (-)) 
converges to dg R (qo, qi). For every n, denote by <p n (-) the flow of X n . Up to the extraction of some 
subsequence, we have X n (-) -> X(-) for the weak topology in L 2 (0,1; V) C L 2 (0,1; T So+e+1 (TM)). 
Let <£>(•) = ip x (•). Lemma 2 from [5] implies that —>• (p{l)\u strongly in H s as n —> +00 

for every compact subset U of M. Since qo has a compact support, this implies that the horizontal 
curve q(-) = (p(-) ■ qo satisfies q( 1) = q±. But 

A(X(-)) < liminf A(X"(-)) = ^d^ R (q 0 , qi) 2 . 

n—> 00 Z 

□ 
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2.5 Sub-Riemannian geodesics and singular curves on shape spaces 

We assume that S is a shape space in M of order 1, and that V s (M), s = so + £, is equipped 
with a strong right-invariant sub-Riemannian structure induced by the Hilbert space (V. (•,•)) of 
vector fields on M, with continuous inclusion V T s+k (TM) for some k ^ 1. 

The differential structure of the space of horizontal systems and the endpoint map¬ 
ping. Without further assumptions on S (for example, the existence of a local addition), it is not 
known whether H 1 ( 0,1; 5) admits a natural differential structure [35]. However, Condition 4 of 
Definition 3 shows that the space of all horizontal systems with starting point go is in one-to-one 
correspondance with L 2 (0,1; V) through the identification X(-) <-> (<p x (-) ■ qo,X(-)). This endowes 
the space of horizontal systems from go with a smooth manifold structure with a single coordinate 
chart. We will always identify it with L 2 ( 0,1; V) through this correspondance from now on. 

Geodesics. Fix an initial point qo and a final point q\ in S. The endpoint mapping from qo is 

end f 0 (*(•)) =V5 A (1) Wo = R qo ° end(X(•)), 

where end(X(-)) = ip x ( 1) £ V S (M). It is of class C k , as a composition of end(X(-)) and R qo , which 
is of class C k on Im(end) C T> s+k {M). A geodesic on S between the states qo and q\ is a horizontal 
system (g(-),X(-)) joining qo and q\ that is a critical point of the action among all horizontal 
systems with the same endpoints, i.e., a critical point of A(X(-)) = A(q(-),X(-)) restricted to 
end” ({gi}). In other words, for any C 1 -family a i->- X a (-) £ L 2 { 0,1; V) such that <p Xa {l) ■ q 0 = qi 
for every a, with X 0 = X, we have 


d a A ( A a (-))| Q=0 = 0 . 

Remark 12. If (q(-),X(-)) is a geodesic, then ip x is the minimal lift of g(-). Indeed, if X(t) does 
not belong almost everywhere to T~L q = (ker^j)- 1 , it admits a unique split X = A 0 + X 1; with 
X 0 (t) £ ker£ q ( t ) almost everywhere and A 0 (-) ^ 0. But then X s = X — sA" 0 generates the same 
horizontal curve t >—> q(t) and 

d s A(X s ) = ^ [ {X o (t),X o (t))dt^0. 

Jo 

Consequently, a geodesic (g(-),A(-)) is uniquely determined by its trajectory g(-), and we can 
identify geodesics with trajectories of geodesics. We will do so in the rest of the section. 

Remark 13. Any sub-arc of a geodesic is also a geodesic (up to a reparametrization). 

In Section 1, we saw that there can be three kinds of geodesics on V s (M), depending on the 
image of the derivative of the couple (A, end). This is also true here, at least in the general case 
where S has infinite dimensions. If (q(-),X(-)) is a geodesic, that is, a critical point of the action 
with fixed endpoints, then one of two possibilities is satisfied: 

1. There exists (A,pi) eRx T*^S \ {0,0} such that 

d endf 0 (*(•))'+ AdA(A(-)) = 0. (9) 

This is a Lagrange mutiplier characterization, which splits into two subcases: 

• If A ^ 0, we can take A = —1 so that dA(X(-)) = dendq o (A(-))*p!, and say that q(-) is 
a normal geodesic. 
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• If A = 0, p\ f 0 and dend^(X(-))*pi = 0, in which case AT(-) is a singular point of 
end^ 0 . We say that q(j is an abnormal geodesic. 

2. The image of (dA(X(-)), dendg 0 (X(-))) is a proper dense subset of R x T qi S: there is no 
Lagrange multipliers characterization. We say that q(-) is an elusive geodesic. 

Indeed, for a geodesic, the couple of linear maps (d^4(X), dend^ o (X)) can’t be surjective. Then, 
if its image is not a proper dense subset (i.e., if we are not in the second case), it is included in a 
closed hyperplane, and we are in the first case. 

2.5.1 Statement of the results 

In this section, we define the Hamiltonian and the corresponding geodesic flow. 

Canonical symplectic form, symplectic gradient. We denote by oj the canonical weak sym- 
plectic form on T*S , given by the formula 

u(q,p).( 6 q 1 , 6 p 1 ; 6 q 2 , 6 p 2 ) = Sp 2 (Sq 1 ) - Sp^dqf), 

with ( 5qi,Spi ) £ T(^ qp )T*S ~ T q S x T*S in a canonical coordinate system ( q,p ) on T*M. 

A function / : T*S —> M, differentiable at some point ( q,p ) £ T*S , admits a symplectic gradient 
at {q,p) if there exist a vector S7 u f(q,p) £ T^ q ^T*S such that, for every 2 £ T^ q ^T*S, 

df( q , P )( z ) = U](X7 u f(q,p),z). 

In this case, this symplectic gradient V“/(g,p) is unique. 

Such a gradient exists if and only if d p f(q,p) £ T**S can be identified with a vector in T q S 
through the canonical inclusion T q S T**S. In that case, we have, in canonical coordinates, 

V“/(g,p) = (d p f{q,p),-d q f(q,p)). 

The normal Hamiltonian function and geodesic equation. We define the normal Hamil¬ 
tonian of the system H s : T*S —► R by 

H S (q,p) = ^p(K q p) = ^p{€ q K v £*p). 

This is a function of class C k . Moreover, for every (q,p) £ T*S and Sp £ T*S , d p H s (q,p)(dp) = 
5p(K q p ), so that we can identify d p H s (q,p ) £ T**S with K q p £ T q M: V w H s (q,p) is well-defined 
on T*S, and given in canonical coordinates by 

V^H S (q lP ) = ( K q p,~d q (K q p)*p ). 

Note that is of class C k ~ 1 on T*S. Hence, if k ^ 2, it admits a unique local flow, called the 

Hamiltonian geodesic flow. 

Theorem 4. Assume k ^ 1. Then a horizontal curve q(- ) is a geodesic if and only if it is the 
projection of an integral curve (<?(•),p(-)) ofV^H. That is, ( q(t),p(t )) = V w H s (q(t) , p(t)) for 
almost every t in [0, 1] . 

Assume k ^ 2. Then V U H admits a global flow on T*S of class C k ~ 1 , called the Hamiltonian 
geodesic flow. In other words, for any initial point (qo,Po) £ T*S, there exists a unique curve 
t ( q(t),p(t )) defined on all o/M, such that (q(0),p(0)) = (qo,Po) an d, for almost every t, 

{q{t),p(t)) = V“H s (q(t),p(t)). 

This curve is of class C k in time and C fe_1 with respect to (qo,Po). 
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We say that p(-) is the normal covector along the trajectoy. 

Remark 14. One can then prove that any normal geodesic is, locally, a minimizing curve on 
small sub-intervals. This is a straightforward generalization of the proof of the same result for 
finite dimensional sub-riemannian manifolds: use the Hamiltonian flow allows the construction of 
a calibration of normal geodesics [1, 42], 

This theorem will be proved in Section 2.5.2, concurrently with Proposition 6 found in the next 
paragraph. Along the way, wc will also prove the following result. 

Proposition 5. Let tp(-) be a horizontal curve from e on T> S (M). Then q(-) = (p(-) -qo is a normal 
geodesic on S if and only if there exists po £ T* o S such that tp(-) is a normal geodesic on T> S (M) 
with initial normal covector Pq = £, qo Po ■ 

Singular horizontal systems. Define the abnormal Hamiltonian Hq : T*S x V —>■ R by 

HS(q,p,X)=p(Z q X) = Z* q p(X). 

For fixed X, the mapping ( q,p ) >->■ Hy(q,p,X) is of class C k and d p Hy (q,p, X) can be identi¬ 
fied with f q X. This mapping therefore admits a symplectic gradient, denoted (g,p, X) £ 

T^ q p )T*S and given, in canonical coordinates, by the formula 

V“R 0 S (g,p,X) = (ZqX,-d q (ZqX)*p). 

This gradient is of class C k ~ 1 in ( q,p ) and linear in p. We have the following result. 

Proposition 6. A horizontal system ( q(-),X (•)) is singular if and only if there exists p : t i-4 
P(t) £ TWiS \ {0} °f Sobolev class H 1 in time such that, for almost every time t in [0, 1], the 

abnormal Hamiltonian equations 

( 4(t),p{t )) = V u H$(q(t),p(t),X(t)) = (£ q ( t )X(t), -d q (£ q{t) X(t)y P (t)) , 
is satisfied, with p(0) ^ 0 and for almost every time t in [0,1] , £ q / t \P(t) = 0, that is, 

£ q (t){V) C ker pit). 


In this case, dend^- 0 )(Al(-)))*p(l) = 0. 

We say that p(-) is a singular covector along the trajectory. This will be proved in Section 
2.5.2, along with the following result. 

Proposition 7. Let </?(•) = <p x (•) be a horizontal curve from e on V S (M). Then the horizontal 
system ( q(-),X( •)) with q(-) = ip(-) ■ q$ is singular on S if and only if there exists po £ T qo S such 
that <p(-) is a singular curve on V s (M) with initial singular covector Pq = £, qo Po- 

A special case of this result was used in [5] to give examples of singular curves on T> S (M). 

2.5.2 Proof of Theorem 4 and Propositions 5, 6, and 7 

The purpose of this section is to prove, using the momentum map, that normal geodesics and 
singular curves on S are exactly those that can respectively be lifted to normal geodesics and 
singular curves on T> S (M). This is a well-known result in the Riemannian framework of Clebsch 
optimal control [24, 38, 45]. 
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Step 1: Momentum formulation. Recall that the momentum map associated to the group 
action of V S (M) over S is the mapping p s : T*S —> T S (TM)* = T~ S (T*M) given by 

P S {<1,P) =poq = C q P- 

Proposition 8. Consider a solution (q(-),p{-)) to the normal Hamiltonian equation {q(t),p{t)) = 
H s (q(t),p(t)) and let X(-) = Kv(,qt t \P{I) and p s (t) = /r 5 (g(t),p(t)) the momentum along the 
trajectory. Then, for almost every t, 


fi s (t) = ad * xw p{t) 

in the weak sense. In particular, <p x (■) is a normal geodesic on V s (M) with initial covector 
P(0) = Cq(o)P(0) and with momentum p(t) = p s (t). 

The same is true for singular curves: if p(-) is a covector associated to a singular system 
(q(-),X{-)), then ip x is a singular curve in T> S {M) with singular momentum p(t) = ia s (t) and 
initial singular covector P{ 0) = )P(0). 

Proof. Fix Y G r s+1 (TM). Then 

/X s (i)(r) = p(t)(5 9 (^ (t) (r))(^ (t) (x(t)))-5 9 (4 (t) (x(t)))(4 (4) (y))) 

= p(t)(K(x(t)),^(r)]( g (t))) 

= P(t)(£q(t)([X(t), F])) = ad* xw p(t)(Y). 

In the normal case, t p(t) then satisfies the necessary and sufficient condition for X(t) = Kypit) 
to be the control of a Hamiltonian geodesic t i-» ( ip(t),P(t )) on T> S {M) with initial covector 
fj,(0) = P( 0). In the normal case, since X(t) = f*^p(t), (p(-) is also the minizing lift of q starting 
at e. 

In the singular case, 1 1 —> p(t) satisfies the necessary and sufficient condition for X(t) to be the 
control of a singular curve 1 1 -» (p(t), P(t)) on V S (M) with initial covector P(0) = /i(0). □ 

We will also need the converse. 

Proposition 9. Let </?(•) = (p x (•) be a normal geodesic (resp. a singular curve) on T) S (M) with 
momentum p(-). Assume that p(0) = £* 0 po for some po G T* q S, and define the horizontal curve 
q (-) = ip(-) ■ <jo- Then there exists p:iG p{t) G T*^S with p(0) = pg such that ( q(-),p (•)) satisfies 
the normal (resp. abnormal) Hamiltonian equations. Moreover, p. s (q(t) , p(t)) = p(t) for all t. 

Proof. Consider p(-) the solution of the linear Cauchy problem p(0) = po and 

p(t) = -d q (£ qit) X(t))*p(t). 

Let l-i S (t) = p s (q(t) , p(t)) for every t G [0,1]. Then, the same proof as in the previous proposition 
shows that 

p s {t) = ad * x{t) p S {t). 

Then p s (-) = p(-)*p s (0) = tp(-)*p(0), and Proposition 4 implies that <p(-)*/a(0) = p(-). Hence, 
p(t) = p s (t) for every t in [0,1]. 

In the normal case, we need to prove that ( q(t),p(t )) satisfy the normal Hamiltonian equation 
on S. In the singular case, we already have p(t) = —dqH^(q(t),p(t), X(t)), but we must also prove 
that £*p)P(i) = 0 on V. 
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Let us start with the case where ip is a normal geodesic. Then X(t) = Kyg(t) = Kyf( (t yp(t) 
for almost every t, so q(t) = ^ q ^X(t) = K q ^p{t). Then 

-\d q {p{t)K q[t) p{t)) = ~{d q {€ qit) K v £* (t) p(t))*p{t) 

= -(d q (Z qit) )KyQ t) p(t))*p(t) 

= -d q (Z qlt) X(t))*p(t) 

= p(t). 

Hence ( q(t),p(t )) does satisfy the normal Hamiltonian equation on S. 

In the singular case, for any Y £ V, £ q / t \P(t)(Y) = g s (t)(Y) = g(t)(Y) = 0 because p x is a 
singular curve with momentum /i(-), which finishes the proof. □ 

Combining those results, we see that solutions of the normal (resp. abnormal) Hamiltonian 
equations on S are exactly those curves that come from normal geodesics (resp. singular curves) 
on T> S (M) with initial momentum of the form f qo Po- In particular, for k ^ 2, the completeness of 
the the normal geodesic flow on T*T> S {M) implies that is a complete vector field on T*S. 

Step 2: Solutions of the Hamiltonian equations are indeed normal geodesics and 
singular curves. Recall that if (q(-),X(-)) is a normal geodesic (resp. a singular system), then 
d-Apf) = dend^ 0 (X)*pi (resp. 0 = d end^ Q (X)*p\) for some p\ £ T*^S. But the following partial 
converse is trivially true. 

Lemma 3. Let (g(-),X(-)) be a horizontal system with q( 1) = q\. Then: 

• If dA(X(-)) = dend^ 0 (X(-))*pi for some p\ £ T qi S, then q(-) is a geodesic (and, therefore, 
a normal geodesic). 

• If 0 = dend^ 0 (X(-))*pi for some nonzero p\ £ T qi S, then X(-) is a singular point of end qo : 
q(-) is an singular curve. However, it may not be a geodesic. 

Remark 15. On the other hand, (dH(X(-)), dend^ o (W(-))) having proper dense image does not 
give any useful information on X(-), so the elusive geodesic cannot be found this way. 

Now, note that end^ 0 = R qo o end, where end : P 2 (0,1; V) —» V s (M) is the endpoint map from 
e on T> S (M): end(X(-)) = p x {l). Consequently, dend^ o (X)* = dend(X)*dl? 90 ((/3 X (l))*. For 
p £ V s (M) and 5p = Y o p £ T v V s (M), we have 

dRq 0 {p).5tp = f,tp.qY, 


so that d R qo (p)*p = (dR^Yfff.gp. 

We deduce from the previous discussion that (<?(•), X(-)) is a normal geodesic (resp. singular 
curve) between go and gi if and only if p x is a normal geodesic (resp. singular curve) whose 
covector P(-) satisfies 

P(l) = (dp- 1 )*^! 

for some p\ £ T qi S. But P(l) = (dP“ 1 )*/x(l) with p(-) the momentum of p, so this condition 
becomes 

^(1) = Cp-qPi = P S (<H,Pi)- 

Proposition 9 then proves that there exists p(t) £ T q{t )S such that ( q{t),p{t )) satisfy the 

normal (resp. abnormal) Hamiltonian equations, which completes the proof. 
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2.5.3 On inexact matching 

It should be emphasized, again, that other geodesics (the so-called elusive geodesics ) may also 
exist. Moreover, finding every singular curves and axtracting the abnormal geodesics from them 
is a daunting task. However, when performing LDDMM methods and algorithms for inexact 
matching, one aims to minimize over L 2 ( 0,1; V) functionals of the form 

WO) = WO) + ff(<r Y (!)) = WO) + 9 ° endf 0 (*(•)). 

In this case, X(-) is a critical point if and only if 

dA(X) = -dendf 0 (X)*dg(q x (l)). 

The trajectory induced by such a critical point X is therefore automatically a normal geodesic, 
whose covector satisfies p{ 1) = —dg{q x ( 1)). This means that one need consider neither abnormal 
nor elusive geodesics when looking for minimizers of J. 

Consequently, the search for minimizing trajectories can be reduced to the minimization of 

p(t){K qW p(t))dt + g(q( 1)) 

among all solutions of the control system 

4(t) = K q (t)P(t), 

where p(-) is any covector along <?(•) and is L 2 in time. 

Remark 16. In other words, we reduced ourselves to studying another sub-Riemannian struc¬ 
ture SRy T „ 5 = (T*S, g K , K), where q i— K q is the vector bundle morphism T*S —> TS and 

9q(PUP2) =Pi(K q p 2 ). 

This is a weak sub-Riemannian structure, whose normal geodesics coincide with those of SRy. 

This reduction is very useful in practical applications and numerical simulations, since, when S 
is finite dimensional, we obtain a finite dimensional control system, for which many optimization 
methods are available. See [6] for algorithms to minimize such a functionnal in the abstract 
framework of shape spaces in R d . 

2.6 Some examples of geodesic equations on classical shape spaces 

In this section, we assume that M = Let So be the smallest integer such that so > d/2. For 
every integer s ^ so + 1, the group V s (R d ) coincides with the set of diffeomorphisms ip of such 
that tp — Id R d £ H s { R d ,R d ), and is an open subset of the affine Hilbert space Id R d + H s { R d ,R d ) 
[19]. For a vector field X , we will also write X = X 1 e\ + • • • + X d ed = (X 1 , ..., X d ), with (e,) the 
canonical frame of R d . 

The Euclidean inner product of two vectors v and w of is denoted by v • w. The notation 
v T stands for for the linear form w v ■ w. Conversely, for a linear form p £ (R d )*, we denote by 
p T the unique vector v in such that p = v T . 

2.6.1 Spaces of continuous embeddings of compact manifolds and Gaussian kernels. 

We consider the Hilbert space V C H s (R d ,R d ) whose reproducing kernel is the Gaussian kernel 

K(x,y)p = e~ L ^ L p r = e(x - y)p T , 
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. M-* 


where we denoted e - ^ - = e(x) for readability. We let TV be a compact Riemannian manifold, 
and S = C°(N, R d ), a shape space of order 1 (and even 0) in R d . We have T*S = S x C°(N,M. d )*, 
with the elements of C°(N, R d )* being identified to maps S —> (R d )* with coefficients in the space 
A i(N) of Radon measures on N, that is, elements of M.(N) (§) (R d )*. 

Then, for q : N — > R rf continuous and p = (pi,. .. ,Pd ) £ M(N) ig 

{K V C q P)(x) = [ e(x-q(a))dp T (a). 

J N 

Hence, we have K q p £ T q S = C°(N, R d ) given by 

(■ K q p)(a ) = [ e{q(a) - q(a'))dp T (a'). 

Jn 


The normal Hamiltonian is 


H S (q,p) = \ f e(q{a) - q{a'))'y]&pi{a)&pi(a') = )- f e(q(a) - q(a'))[dp{a) ■ dp{a')}. 
z Jn 2 1 Jn 2 

Using e(x — y) = e(y — x) and de(x).v = —j;(x ■ u)e(x), we get for any Sq £ T q S , 

d q H s (q,p).Sq = -- f [( q{a ) - q(a')) ■ Sq(a)]e(q(a) - q(a'))[dp(a) ■ dp(a')]. 
f 7 J N 2 

In the end, the geodesic equations read 

q(t,a)= / e(q(t,a) - q(t,a'))dp T (t,a r ), 

Jn 

p(t,a) = - f e(q(t,a) - q{t,a')){q(a) - q(a')) T [dp(t,a') -dp(t,a)}. 
a Jn 

In the special case of landmarks, we can write S = (R d ) ra , writing q = (xi,... ,x n ) and p = 
(pi,... ,p n ) £ (R d *) n , we get (omitting time for readability) 


= Y1 


e 2 ” p, , 


j=i 




i = 1 ,..., n. 


Pi = y2 e 2 (jpi ■ Pj){xi - Xj) 1 


j =i 

Remark 17. It is worthwile to note that in the neighbourhood of q = (x±, ..., x n ) such that i ^ j 
implies Xi ^ Xj , the structure here is actually Riemannian instead of simply sub-Riemannian. 

2.6.2 An action of higher order 

In this paragraph, we keep our Gaussian kernel K as above and N our compact Riemannian 
manifold, but we take the shape space S = TC°(N , R d ) = C°(N , R d ) x C°(N, R d ), a shape space of 
order 1 for the action 


ip ■ (q,v)(a ) = (ip(q(a)),dip(q(a)).v(a)), a £ N. 
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The three conditions for a shape space of order 1 are easily seen to be true. The infinitesimal 
action of a vector field X on is 

Zq,v{X){a) = (X(q(a)),dX(q(a)).v(a)), a G N. 

Now a momentum on S is a couple (p, l) of maps N — > (R“)* with coefficients in M(N), and 

tq,v(P> l )( X ) = [ M°-)(X(q(a)) + f dl(a)(dX(q(a)).v(a)). 

Jn Jn 

Therefore, 

rV C q ,v(Pi 00) = [ e i x ~ l{a))dp T (a) e(x- g(a))[v(a) • {x - q(a))]dl T (a). 

Jn a Jn 

Writing £ qtV K v £* tv (p, l) = (w 1 ,w 2 ), and e a , a > = e(q(a) - q(a')), we get 



The normal Hamiltonian is given by 

H S (q,v,P,l) = \ [[ e a , a -[dp(a) ■ dp(a')] - T / e a y[v(a') ■ v{a)][dl(a) ■ dZ(a')] 
z J Jn 2 Za Jn 2 

~ 7 ^ JJ 2 e a ,a'[(v(a)) • {q{a) - q(a'))}[dl(a) • dp(a') + dp{a) • d l(a')] 

+ ^ JJ N2 e a,a’[v(a') • ( q(a ) - q{a'))][v{a) • ( q(a ) - q(a'))][dl(a) • d l(a')]. 

From there, the geodesic equations are easily deduced. For example 

i(a) T =-d v H s (q,v,p,l) T (a) =- [ e a , a /v(a')[l(a) • dl(a')] 

a Jn 

+ 7 T ( e a ,a'(<?(a) - q(a'))[l{a) • d p(a') + p(a) • dZ(a')] 

Jn 

- ~2 [ e a ,a>(q(a) - q{a'))[v(a') • ( q(a ) - q(a'))][l(a) • dZ(a')]- 
a Jn 

In the case of landmarks, we take q = (x±, ..., x n ), v = (vi ,..., v n ), p = (jpi ,... ,p n ), and 
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I = (7i,..., l n ). Denoting Xij = Xi — Xj ,we get 

X I [Pj - ^i v 3 ■ e ( x ij)> 

3=1 ' ' 

n ^ ^ 1 1 \ 

vf ( -~l v i • x i,j]Pj ~ ~2\ v i ■ ( x i,j)][ v j ■ x ij] l J + ~i v i ' v j] l j ) e ( x i,j), 

3 =1 ' ' 


h — n ^ ( l^i ■ y e ( a 'j,i)' i; .7 “I" 2^. ^ ( [^» ' Pj liPj\ e { x i,j) x i 


3 =1 


i=i 


J=1 

The derivative of is slightly more complex, as 
1 


Pi 




j=i 


Pi ' Pj ['Ci ' S-ij] |Pj ‘ ‘ ^j] [^i ' Pj][^i ' ) e { x i,j) x i,J 


1 " 

— ^2 ^ I ([p® ' x i,j][li ' Pj ] — [P? ' X i,j][Pi ■ y ) e ( x i,j) x i,3 

3 =1 

1 " 

— ~3 (|p® ■ a '®d]['P? ■ ■ *j]) e ( x i,j) x i,j 


3=1 


+-Y 1 ■ pj] w ® - ^ • y w j) 


i = 1 ,... ,n 


i = 1,..., n. 


i=i 

Remark 18. Again, in the case of landmarks, there is an open and dense subset of S on which 
the sub-Riemannian structure is actually Riemannian. 


2.6.3 The case of sub-Riemannian Gaussian kernels 


Here, we compute the sub-Riemannian Hamiltonian geodesic equations for the shape space S = 
C°(N, R d ), with N a compact Riemannian manifold, and with V defined by the reproducing kernel 

i*- h i 2 k k 

K(x,y)p = e ^2p(X r (y))Xr(x) = e{x - y)'Y^ J p{X r (y))X r (x), 

r=l r—1 

for Xi, ..., Xk smooth bounded vector fields on R d with bounded derivatives at every order, linearly 
independent at every point x of R d . Note that any X G V is everywhere tangent to the distribution 
of subspaces of dimension k generated by the X r s. 

Now for q G S and p = (pi,... ,p d ) G T*S = M{N) ® (R d )*, we get 

{K v £*p)(x) = V [ e(x-q(a))dp(a)(X r (q(a)))X r (x), 

r—l 


where dp(a)(X r (q(a))) G M(N) is defined by X^(q(a))dpi(a) + ■ ■ ■ + Xf(q(a))dpd{a). 
Then we get 


K q p(a) = K v £*p(q(a)) 



e{q{a) - q(a'))dp(a')(X r (q(a')))X r (q(a)), 
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and compute the reduced Hamiltonian 

k 


H S (q,p) = ^ E _// 2 e{q(a) - q(a , ))dp(a)(X r (q(a)))dp(a')(X r (q(a'))). 

Then, reintroducing the notations e a , a ' = e{q(a) — q{a')), we obtain 

k 

d q H s (q,p).Sq = -- V [[ e ata ’{5q(a) ■ ( q(a ) - q{a'))}dp(a)(X r (q(a)))dp(a')(X r (q(a'))) 
k 

E/7 e a ,a'dp(a)(dX r (q(a)).5q(a)dp(a'){X r (q(a'))). 

, J Jn 2 


r— 1 ' 

In other words, as a Radon measure, 

k f. 


•\\T 


| %H s (q,p){a) = E / e a,a , dp(a / )(X T .(g(a / ))) ( p(a)(dl r (g(a))- p{a){X r {q{a))){q{a) - q(a')) 

f-[JN V °- 


Here, p(dX r (q)) in A4(N) ® (R d )* is the (R rf )*-valued Radon measure given by 

dp(a)(dX r (q(a)) = dpi (a) ®dX*(g(a)) H-hdp d (a) ® dX^(g(a)). 


Therefore, the Hamiltonian geodesic equations are given by 
k „ 

9(a)=E / e a,a>dp(a')(X r (q(a')))X r (q(a)), 


r —1 
fc 


'TV 


P(a) ='£j N e a ,a'dp(a')(X r (q(a’))) ^p(a)(X r (q(a))){q(a) - q(a')) T - p(a)(dX r (q(a)))^ . 

In the case of landmarks, q = {x \,..., x n ) £ (M d ) ra and p = (p lt ... ,p n ) £ (IR d ) n *, and denoting 


Xi — Xj = , we obtain 


k n 

Xi = ^ ^^ c{xjj^pj( < X r {xj'y^X r {xi)^ 
r= 1 j=l 

k n / ^ 

Pi = EE e(xij)pj(X r (xj)) ( pi(X r (xi')')Xij pi(d^f r (x^)) j . 

r _i j—i V 


( 10 ) 


Remark 19. An interesting fact is that e(xij) —> Sij as a —>• 0, where 5jj = 1 if i = j and 0 
otherwise. Hence ^ e(xij)xij —> 0 as a —> 0 because Xu = 0. Therefore, as a goes to 0, the geodesic 
equations become 

k 

Xi — ^ ) Pi (,X r (Xj ) ) X r (Xj ), 

r =1 
k 

Pi — ^ ) Pi (X r {Xi Y)Pi (dX r (a^i)). 

r= 1 

For each couple t n- (xi(t) , Pi{t)) £ T*R d , we recognize the hamiltonian geodesic equation for the 
sub-Riemannian structure A induced on R d by Xi,..., Xk [42] : the dynamic is that of a system 
of n points without interaction in (IR d , A). 

So Equations (10) can be seen as a perturbation of such a dynamical by creating an interaction 
between the points Xi, which decreases exponentially with respect to the square of the Eucliddean 
distances between the points. 
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3 Equivariant mappings and lifted shape spaces 

The general purpose of this section is to define and study the applications of equivariant mappings 
between shape spaces. In particular, we will use them to study symmetries in both Hamiltonian 
flows, and to consider ILH shape spaces of infinite order. Another use for equivariant maps is the 
introduction of lifted shape spaces, which allow the consideration of additional information along 
deformations, and are always strictly sub-Riemannian even for finite dimensional shape spaces. 

3.1 Equivariant maps 

We keep the notations from the previous sections. 

Definition 6 . A mapping tt : S > S between two shape spaces of respective order £ ^ £ is said to 
be equivariant if for every diffeomorphism p and state q £ S, <f> ■ tt( q) = tt (tp ■ q). 

The composition rule for derivatives shows that for every X £ T So+e (TM), we have 

dlT^qf^qX f,7r(q)X- 

In particular, a horizontal curve on S is projected to a horizontal curve on S with the control X(-). 
Let us give a few examples. 

Example 5. Additional landmarks: Take n ^ n £ N*, and S = Lmk n (M) and S = Lmkn(M). 
We define tt : Lmkh(M) Lmk n (M) by 

tt(xi,. ..,x a ) = (xx, .. .,x n ). 

The mapping tt is an equivariant submersion for the diagonal action of T> So+1 (M): for every tp in 

V So+ 1 (M), 

Tr(ip ■ (xi, . . . , Xn)) = Tr(ip(xi), . . . , lf(Xn)) = . . . , p{x n )) = P ' TT^Xl, . . . , Xn). 

This is the mapping used in [54] for the diffeon method. 

Example 6 . Embeddings of submanifolds: Let A be a smooth compact Riemannian manifold, 
N C N a submanifold of N without boundary, and let S = and <S = C°(S,M). Then 

the map 7 r : < 71 —>• < 7 |tv is an equivariant submersion. 

This was used in [ 6 ] to prove that the LDDMM methods on discrete shapes converge to the 
continuous case as the discretization gets finer and finer. 

Example 7. Action on the tangent space: Let S be a shape space in M of order l. Then 
S = TS is a shape space of order £ + 1 for the action of T> So+e+1 (M) given by the differential of 
the action on S 

<P ' (<L v) = (p- q, d q (ip ■ q)(v)) £ T v . q S. 

Then the projection tt( q, v) = q is an equivariant submersion for these actions. 

In the case of landmarks, i.e. S = Lmk n {M ), we can write (q, v ) = (aq,..., x n , iq,..., v n ), 
with Vi £ T Xi M , and 

P- {q,v) = (p(x 1 ),...,p{x n ),dp Xl (v 1 ),...,dp Xn {v n )). 

Example 8 . Inclusion in a shape space of lower order: Consider 5 = Emb a ° +e {N,M) and 
S = Emb“ 0 + ^ + 1 (A, M), with N a smooth compact manifold. Then the inclusion S ^ S is trivially 
equivariant. 

Depending on the properties of 7 r, many things can be deduced from the existence of such a 
mapping. 
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3.2 Equivariant mappings and the geodesic flow 

We keep the notations and setting of the previous section, and take n of class C°°. Fix V a Hilbert 
space of vector fields with continuous inclusion in T So+e+k (TM) for k ^ 1. Then a horizontal 
geodesic on S can be pulled back to a geodesic on S through n. 

Proposition 10. Let (g(-),p(-)) be a normal geodesic (resp. a singular curve) and the corre¬ 
sponding covector on T*S, such that q{ 0) = 7r(<jo). Let A(-) be the corresponding vector field in 
L 2 (0, 1;E). Then for q(-) = <p x {-) ■ % and p(-)) = dn(q(-)*p(-), (q(-),p(-)) is a normal geodesic 
(resp. a singular curve). 

Conversely, let (<?(•), p(-)) be a normal geodesic (resp. a singular curve) on S, and denote 
q(-) = 7r (<?(•))• If p(0) = dn(q(0))*po for some po G Ti g S, then q(-) is a geodesic (resp. an 
abnormal curve) with initial covector po. 

Proof. Let ( q(-),p( •)) be a normal geodesic (resp. a singular curve) and the corresponding covector 
on T*S , such that q( 0) G n(q 0 ). Define ( q{-),p{ •)) as in the hypothesis. We just need to prove that 
p, s (q(t) , p(t)) = p s (q(t) , p(t)) thanks to Propositions 8 and 9. But since dTr(q).fq = ^(q) we have 
C(g) = so for all t G [0,1], we do get 

H S (q(t),p(t)) = Q t) p(t) = f* m p(t) = p s (q{t),p(t)). 

Conversely, let (<?(•),p(-)) be a normal geodesic (resp. a singular curve) on S such that p(0) = 
dTr(q(0))*p o for some p 0 G T| o <S, and denote q(-) = n(q(-)). Then the flow of the control X(-) 
corresponding to q(-) is a geodesic (resp. singular curve) from e on V s (M) with initial momentum 
^|(o)^(b) = and therefore induces a normal geodesic (resp. singular curve) on S with initial 

covector po- But since for every t in [0,1] 

p x (t) ■ q( 0) = n(ip x (t) ■ q(0)) = n(q(t)) = q(t), 


we proved the converse. □ 

This result can be used to find symmetries in the geodesic flow of S, for example when S is the 
quotient of S with a group action that commutes with that of V So+l (M). Indeed, it shows that 
the covector p(-) stays of the form dir(q(t))*p(t) for all t as long is had this form at t = 0. 

When 7T is an immersion, can also use it to find elusive geodesics on S , and also consider ILH 
shape spaces of infinite order. 

3.2.1 ILH shape spaces of infinite order 

An ILH shape space of infinite order is an Inverse Limit Hilbert space S°° (see [35, 43]), inverse 
limit of a decreasing sequence <S* +1 C S\ with each S l a shape space of increasing order li > £j_i, 
such that the inclusions <S* +1 S l are continuous, equivariant and have dense image. As a 
consequence, acts smoothly on S°°. 

Example 9. A typical example is S°° = Emb°°(S', M), the manifold of all smooth embeddings 
of a compact manifold S in M. It is the inverse limit of the sequence of Hilbert manifolds 
Emb“ 0+l (S', M), iGN with a 0 > dim(5)/2. 

When V is a Hilbert space of vector fields with continuous inclusion in every T S (TM), it defines 
a sub-Riemannian structure of class C°° on each 5®, and also on S°°, through the infinitesimal 
action. 
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There is no problem with defining the normal and abnormal Hamiltonians on 1 S 00 as usual, but 
the proofs of Section 2 do not work: linear differential equations may not have a unique solution 
on ILH spaces. However, it is very easy to check that, for fixed qo £ <S°° and p £ T* Q S l , we also 
have p £ T* Q S l+ i for every j ^ 0 and 

H S \q,p) = H si+j (q,p), V" H S \q,p ) = V“ H 5 ^ (q,p). 


Consequently, H s °° (q,p) = H s \q,p ) and H s °° (q,p) = H s \q,p). Therefore, (g(-)iP(')) is a solu¬ 
tion of the normal geodesic equation in T*S°° if and only if it is a solution to the normal geodesic 
equation in T*S l , and we know that this last equation has a unique solution defined on all of K.. 
Since T* o S°° = L>i^oT* a S' 1 . we obtain the existence of a global geodesic flow on T*S°°. 

This discussion shows that, on an ILH shape space of infinite order, the Hamiltonian geodesic 
flow is still complete, and we can solve the geodesic equations as long as we can solve them in each 
S\ 


3.2.2 Finding symmetries in the normal geodesic equations 

In this section, we give a method to investigate symmetries on a shape space S on which a group 
G acts on the right with an equivariant action with respect to that of V s (M). The goal is to prove 
that the momentum for the action of G is constant along normal geodesics and singular curves, 
in a way that is similar to what happens with classical Clebsch optimal control problems with 
symmetries [24, 38, 45]. The problem is that in our setting, we need to consider actions that do 
not come from Lie groups: G will usually itself be a group of diffeomorphisms of fixed regularity. 
We will solve this problem in the most important example in a way that can easily be generalized 
to a wider range of problems. 

We will only consider shape spaces of the form S = Emb a ° +e (N, M), with t ^ 1, TV a smooth 
compact manifold, and ao the smallest integer greater than dim(iV)/2. The action of T> S (M) is 
given by the composition on the left: p ■ q = p o q 1 so that f q X = X o q. We want to consider the 
symmetry group G = T> ao+e (N), which acts on S contnuously on the right by right composition: 
q ■ r/ = qo i) for q in S and p in G. This action is indeed equivariant: p o (q o p) = {p o q) o p. Note 
that G is a Hilbert manifold but only a continuous group, and its action S is only continuous. 

However, if q £ Emb“ 0+ ^ +1 (iV, M), the omega lemma (see [26, 43, 47]) implies that p n- q o p 
is of class C 1 . Consequently, we have an infinitesimal action 

C : T e G = T ao+ \TN) T q S 

given by <^ q (Y)(a) = d q(a).Y(a) for every a £ S. This allows us to define a momentum map at q 
for the action of G by 

m(q,p) = C q P = d q*P e T* e G = Y- a °-\T*N) 

Fix V a Hilbert space of vector fields on M that has continuous inclusion in T So+e+k (TM) for 
k ^ 1. Also fix an initial point qo £ S, and assume that qo £ Emb“ 0+£+1 (A r , M). Let (q{-),p(-)) 
be a horizontal geodesic (resp. an abnormal curve) for the sub-Riemannian structure induced by 
V on S, and X(-) the corresponding control. Then, for every time t, q(t) = p x (t) o qo- But since 
X £ L 2 ( 0,1; V), p x (t) is of class H So+e+1 , so q(t) £ Emb“ 0+f+1 (./V, M). Therefore, the momentum 
m(t) = m{q(t),p(t)) with respect to the action of G is well-defined along the trajectory. 

Proposition 11. Let ( q(‘),p (•)) a horizontal geodesic or an abnormal curve for the sub- 

Riemannian structure induced by V on S with covector p(-), and A(-) the corresponding con¬ 
trol. Assume that q(0) belongs to the dense subset Emb“ 0+f+1 (IV, M). Then the momentum 
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m (t) = m (q(t),p(t)) with respect to the action of G is well-defined and constant along the trajec¬ 
tory: m(t) = m(0). 

In particular, if k ^ 2, this momentum is constant under the Hamiltonian geodesic flow ofT*S 
as long as the starting point go belongs to the dense subset Emb a ° +e+1 (N, M). 

Proof. Fix Y G T ao+e {TN) = T e G. Then 

m (t)(Y) = p{t)(dq(t).Y). 

Then m(-)(Y) : [0,1] —> R is of class H 1 , because p(-) and dq(-) are of class H 1 . This is trivial for 
p{-), and true for dq(-) because (g(-), dg(-)) is the deformation induced by X(-) in the shape space 
TS of order £ + 1 of (qo,dqo). Moreover, we have d t dq = dX(t)(q(t)).dq(t). Therefore, for almost 
every t in [0,1], 

d t m(t)(Y) = p(t)(dq(t).Y) + p(t)(dX(t)(q(t)).dq(t).Y). 

But in canonical coordinates, since £ q ( t )X(t) = X(t) o q(t), 

p(t)(dq(t).Y) = -d g (X(t)oq(t))*p(t)(dq(t).Y) = -p(t)(d q (X(t)oq(t)).dq(t).Y) = -p(t)(dX(t)(q(t)).dq(t).Y) 

and we do get d t (m(t)(Y)) =0. □ 

Remark 20. Note that when considering S = Emb°°(]V, M), the action on the right by V°°(N) is 
smooth and has a well-defined momentum map. Then the results of the previous section immedi¬ 
ately combine with Proposition 11 to prove that this momentum is preserved under the Hamiltonian 
geodesic flow. 

Remark 21. The general case for abstract shape spaces would be as follows: assume a topological 
group G with a Banach manifold structure acts continuously and equivariantly on the right on a 
shape space S. Assume that go is such that p G G g Q • p is of class C 1 , so that we can define 
a momentum map m at go- Then, along any normal geodesic or singular curve q(j from go with 
covector p(-), the momentum m(q(t),p(t)) is well-defined and constant. 

Application: inexact matching between shapes up to reparametrization. The main 
application of Proposition 11 is to study problems that do not depend on a reparametrization qop 
of the state g. For example, assume we want to minimize 

J(A(.)) = A(A(.))+g(g(l)) 

over all horizontal systems (q,X), where g : S —> R is of class C 1 and g(q) only depends on 
Im(g) (for example, g could be a norm on currents [28, 21]). This is equivalent to saying that 
g(q ■ p) = g(q) for every reparametrization p: if Im(g)=Im(g / ), then g = q' o (g ,_1 o g). We know 
that the trajectory g(-) of every critical point of J is a normal geodesic whose covector p(-) satisfies 
p{l) = —dg(q(l)). 

Now, if go is in Emb a ° +i+1 (N,M), then m(g(l),p(l)) = m(g(l), — dg(g(l))). But since g is 
invariant under the action of G, we have m(g(l), — dg(q(l))) = 0. Therefore, m(g(t),p(f)) = 0 for 
every t, since the momentum is conserved along a normal geodesic. 

Hence, when looking for minimizers of J, it is enough to consider the control system 

q{t) = K q ( t) u{t ), 

where u(t) G T*^<S and m(g(i), u(t)) = 0 for almost every t. The cost becomes 

u (')) = \ J o u(t)(K q{t] u{t))dt + g(g(l)). 
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Morover, m (q(t),u(t)) = 0 if and only if u(t), identified with a vector field (with distributional 
coefficients) along Im(( 7 (f)) thanks to the Riemannian metric on M, is everywhere orthogonal to 
the submanifold We recover the well-known (but never proved in this general setting) 

and widely used fact that, when matching submanifolds, the covector must be orthogonal to the 
surface [13, 39] along a minimizing deformation. 

Remark 22. Exact matching of submanifold is much more complex, because we may again see the 
appearance of abnormal and elusive geodesics. However, we can still obtain a global Hamiltonian 
geodesic flow from any initial point qo £ Emb a ° +f+1 (iV, M), simply by taking an initial covector po 
with m(qo)Po) = 0. This will yield curves that are critical points of the action among all horizontal 
curve with endpoints in the same G-orbit. 

3.3 Equivariant submersions and lifted structures 

The structures defined in Section 2 are usually not Riemannian for an infinite dimensional shape 
space S. For example, for S = Emb“ 0+f (N, M), we require vector fields that are more regular than 
the embeddings we are studying so we cannot obtain the whole tangent bundle of S. However, in 
most practical applications, the shape spaces are discretized into landmark spaces, for which the 
sub-Riemannian structure induced by a space V of vector fields is, itself, Riemannian. 

Sub-Riemannian geometry still has many practical uses in shape analysis. In [6], shape spaces 
with constraints that are linear with respect to the control were studied, which is precisely the 
framework of sub-Riemannian geometry. On the other hand, in [54], Younes gave a new kind 
of LDDMM methods using what he called diffeons, in which the momenta used to induce the 
vector fields deforming the shape have support outside of the shape, but those supports are still 
transported by the deformation. This leads to a sub-Riemannian setting as well. In this section, we 
will establish a framework encompassing the diffeon method in a broader setting using equivariant 
submersions. 

The main use of this setting is to keep track of additional information on a given shape space, 
or to emphasize the importance of certain points in the shape. The resulting normal geodesic 
equations are more complex. We give some examples at the end of the section. 

Let us start with a simple case to clarify the concepts we wish to introduce. Let S and S' 
be two shape spaces of respective orders t and £'. Then S = S x S' is a shape space of order 
1 = max(f, £') for the diagonal action of D So+ \M) 

■ (qiq') = (v? ■ q,<p ■ q')- 

Fix fc > 1 and a Hilbert space of vector fields V with continuous inclusion in T So+e+k (TM). V 
induces on each of S , S' and S a sub-Riemannian structure SRy , SRy , SRy , as described in 
Section 2. We will denote by £ all infinitesimal actions on these three spaces, the index making it 
clear which action is being used. We can then study simultaneous deformations of qo in S and q' 0 
in S' by considering the product space S. 

A different interesting problem to consider is the following. Fix ( qo,q ' 0 ) £ S x S ', and consider 
a horizontal deformation q(-) of qo in S. What kind of deformation of q' 0 does q(-) induce? 

More precisely, q(-) can be lifted to a unique minimal curve <p x in V s (M) with minimal action. 
This lift is the flow of the vector field X(t) = ^~^q(t), where : £ q (t)(V) — > (ker^( t ))- L is the 
pseudo-inverse of £ g ( t ). It induces, in turn, a horizontal deformation 

q'(t) = {q(t), q\t)) = (ip(t) o q 0 , ip(t ) o q ' 0 ) 

of q'o in S ', with 
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This yields a horizontal curve g(-) = in S with 

fa) = tm^ t) fa)- 

Remark 23. The curve g(-) is, in fact, the only one from go whose action in the sense of Remark 
11 coincides with that of g(-). In other words, it is the unique minimizing lift of g(-) to S through 
7 r from go- The converse is trivially true: if 7r(g(-)) = g(-) and their action in the sense of Remark 
11 coincide, then q(t) = almost everywhere. 

We obtain a new sub-Riemannian structure SRf-^ = (£,g 5 ,x), with 

£(q,q') = 

9(q,q')( U ’ V ) = (tq lu ’£q lv )’ 

and 

X(q,q') u = £(g ,q')£q U - 

Remark 24. This structure is particularly hard to study in general, because £“1 is very hard to 
compute, and may not be even continuous in q. This is why our results will assume that £ q {V) is 
closed in T q S for every g, which will mean that £ q (V) = K q (T*S). 

This is a particular case of a lifted sub-riemannian structures on shape spaces through an 
equivariant submersion. We will see the general case in the next section. 

We keep the notations and setting of the previous section. Going back to the end of Section 
3.1, we could define a lifted sub-Riemannian structure on S by lifting that of S through n. The 
vector bundle part would be 

£q = £ir(ij)(^0) 

the vector bundle morphism at q would be given by f, q ^~^ q y and the squared norm of u G £ q would 

be (£„, ^(g)" 1 ^). However, this structure is not easy to study since (;~ l is hard to compute and 
may not even be continuous in g. 

However, in the particular case where £ q (V) is closed, Lemma 2 implies that £ q (V) =Im (K q ), 
where K q = £, q Kyf;*. This allows us to define the following lifted structure. 

Definition 7. Lifted sub-Riemannian structure: Let S and S be shape spaces in M of re¬ 
spective orders £ ^ £. Fix k ^ 1, V a Hilbert space of vector fields with continuous inclusion in 
r so+ \TM), and 7 r : S —>■ S an equivariant map of class C°°. We assume that for every q £ S, 
£ q {V) is closed in T q S. 

The lifted sub-Riemannian structure on S is SR W = (tt*TS, g K , x), where: 
f (t r*TS) q = T; (4) S 
< gf{ui,u 2 ) = ui{K^ q) U2), ui,U 2 £ 

{ Xq=ZqKv£ (q y 

We denote A^ the corresponding action: 

1 f 1 

A n(q(-),u(-)) = - u{t)K n(m) u(t))dt. 

Remark 25. The assumption that £ q (V) be closed usually (although not always) implies that S 
is finite dimensional, for example a space of landmarks. 
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Remark 26. We denote elements of T*S by the letter u instead of p in order to stress that u(-) 
does not satisfy any of the Hamiltonian equations given in the previous section. The covector «(■) 
is simply a control, that can take any value in T*S. 

This structure is equivalent to the one we were looking for, as shown in the following proposition. 

Proposition 12. A curve q(-) £ H 1 (0,1;<S) is a horizontal for SR W if and only if q(-) = <p x (-) ■ 
n(q(0)), with X(t) = ?Wg(t))9(0- Moreover, its action is equal to that of the minimal lift to 
V S °+\ M ) of n(q(-)). 

Proof. Denote q(-) = 7r(g(-)), qo = g(0) and qo = q( 0). Now, assume q(-) is horizontal for SR^, so 
that, there exists u : t TW<S in L 2 such that for almost every t £ [0,1], 


q{t) = ^q{t) K vC q (t)U{t). 

But then 

q{t) = dTr{q(t))4 m K v Q t} u(t). 

Since dn(q(t)).^ t ) = £„($), we get 

<i(t) = iq(t)K V C q (t)U{t) = K q{t) u{t). 

But then Lemma 2 states that if X(t) = Ky^*^u(t), then X(t) = so we do get that 

q(t) = £q£qmq(t) and q(t) = ip(t) ■ qo- Moreover, using the discussion from Section 2.3, we get 

= \ J q u{t){K q(f) u{t))dt = «(■))• 

For the converse, assume that for X(-) = <?(•), we have ip x (•) • qo- Since for every t in [0,1], by 

hypothesis, £g(t)(V) is closed, we know that there exists u(t) £ TW<S such that X(t) = Kvf*^u(f). 
Moreover, 

[ u(t)(K q ^u(t))dt = [ (X(t),X(t))dt. 

Jo Jo 

We just need to prove that q(t) = f q (t) x vf*^u(t) for almost every t. But 


im K vCq(t) u( f) = £, q (t) x (t) = q{t)- 

□ 


Let us give a few examples. 

Example 10. Take n ^ n £ N*, and S = Lmk n (M ) and S = Lmkh(M). We define 7r : 
Lmkfi,(M) —»■ Lmk n (M), as above, by ... ,xa) = (xi,.. ., x n ). Now let V be a hilbert 
subspace of vector fields of class at least H So+1 , with associated kernel K{x,y) : T*M —> T X M. 
Let q = (xi, ..., Xff) and q = n(q). Now take some 

p = (p!,... ,p n ) £ T*Lmk n (M) = T* X M x • • • x T^M. 


Then for a; £ M, 


X(x) = K v C q p(x) = Y^ K {x,Xi)pi, 


31 


so we get 

( n n \ 

y^if(xi ,xj)pj,... ,y^K(xfj,xi)pi I. 

2=1 2=1 / 

On the other hand, we have p = chr*p = (pi,... ,p„, 0, ... , 0), so we also get 

/ n n \ 

KqP = ^2K(X!,Xi)pi,. . . ,^2K(Xn,Xi)pi . 


Example 11. Take n £ N*, S = Lmk n (M) and S = TLmkh(M). We take ir : TLmkn(M) —>• 
Lmk n {M ) the usual projection. Now let V be a hilbert subspace of vector fields of class at least 
H So+ 2 , with associated kernel K(x,y) : T*M —> T X M. 

Let q = (q, v ) = (x±,..., xa, Vi ,..., u„), where <7 = (* 1 ,..., x n ) and Vi £ T Xi M. Now take 
some p = (pi,...,p n ) £ T*Lmk n (M). Again, we have p = d-R*p = (pi,... ,p n , 0,..., 0) £ 
T*TLmk n (M), so we obtain 

Kgp= (w 1 ,...,w n ,w[,...,w' n ), 

with 


n 

Wi = y ^K(xj,Xj)pj £ 

1=1 



which is well-defined as the derivative at Xi of the map K v £*p : M —>■ TM applied to Vi. 

Lifted shapes can be used for studying how the deformation of a given shape induces the 
deformation of a ’’bigger” shape, i.e. one with more information. It can also be used to induce a 
different kind of deformations, such as a deformation through diffeons [54]. 

We can then define horizontal curves for this new structure, the action A and length L of such 
curves (which coincide with the action and length for 7i s ), and the corresponding sub-Ricmannian 
distance df on <S. 

We have the following result, a straightforward consequence of Proposition 12. 

Proposition 13. Let q : [0,1] —> S be a curve of class H 1 , and define <7 : t H> n(q(t)) its projection 
onto S. Then 

1. If q is a geodesic (resp. a minimizing cuvre) on S, so is q, with respect to both SRy and 
SR„. 

2. As a consequence, tt is a lipshitz-1 map: for every < 71,(72 £ M, <^(< 71 ,( 72 ) dg R (n(qi),Tr(q 2 )), 
where df is the distance associated to SR^. 


3.4 The lifted Hamiltonian geodesic equations 

In this section, we compute the normal Hamiltonian for a lifted structure, which will give us the 
normal geodesic equations. We keep the notations and setting of the previous section. Moreover, 
we make the following assumption: 

(Al): The structure SRf- is Riemannian, that is, £ q (V) = T q S for every q in S. 
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This assumption implies that K q is a linear isomorphism. In most practical applications, this 
restricts us to the case where S is a finite dimensional manifold such as a landmark space, but this 
was already a consequence of assuming £ g (V) closed in T q S. 

The Hamiltonian H v : T*S —> S is given by 

( 1 1 

H n (q,p) = max p(£qK v Z* { q)u) - -u^K^u) . 

lit J -n O 

■*{q) y 

V =H^(q,p,u) ) 

Since, for fixed (q,p) € T*S, the right-hand side is strictly convex in u thanks to (Al), the 
maximum is reached if and only if d u Hi(q,p,u) = 0, that is, if and only if 


^7r(g) £,q P A-7r(q)^' 


But, because of Assumption (Al), this is equivalent to 

u = u(p , q) = K~^ (q) K v CqP. 

We obtain for the normal Hamiltonian 

H w {q,p) = H?(p,q,u(p,q)) = ^u{p,q)(K n(q) u{p,q)). 

Now computing the symplectic gradient of H n seems complicated, but when we use the fact that 
d u Hi (p,q,u(j>,q)) = 0 it gets much simpler. Considering that AA(q) = dn(q)Kqd,Tr(q)* 3 , we get 


f d p H 7r (q,p)=^K v ^* ( ~ ) u(q,p), 


1 -d q H*(q,p) = —(dqKqdir(q)*u)* 

p- i dn(q)*u 


+ dq[dTr(q)*u]*K^p. 


Note that we took u = u(q,p) for short and that none of its derivative are involved. This symplectic 
gradient is of class C fc_1 and therefore, if k ^ 2, possesses a well-defined flow. 

Remark 27. We did not prove that integral curves of this symplectic gradient project to sub- 
Riemannian geodesics. This is obviously true is S is finite dimensional, since this is a well-known 
result in sub-Riemannian geometry. The case where S is infinite dimensional is harder, because 
there is no natural differential structure on the space of horizontal systems if S does not possess a 
local addition [35]. 

3.5 Some examples of lifted geodesic equations 

We consider some simple cases of lifted shape spaces and compute their geodesic equations. Note 
that, contrarily to examples from Section 2.6, even the landmark cases will be strictly sub- 
Riemannian. 

Extra landmarks with the Gaussian kernel Let n,m £ N*, and consider S = Lmk n { R d ) 
and S = Lmk n+m (R d ), with n(xi,... ,x n+m ) = (aq,..., x n ). Then T*S = Lmk n ( R d ) x ( R d ) n *, 
and T*S = Lmk m+n (R d ) x (R d )( m+n )*. 

3 Recall that 
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We assume that the reproducing kernel K of V is the diagonal Gaussian kernel, as introduced 
in Section 2.6. Our goal is to compute the geodesic equations for this particular example. In other 
words, we first need to solve the linear system 

K^)U = C^)K v CgP- 

Then, we shall compute dpH^(q,p, u) = £qK v and dqH™ (q,p,u) = (p—^dn^u)((dqKq)dn^u). 

Note that the term (p — dTr'~u)(Kq((dqdTr^)u) does not appear in dqH^(q,p,u) because dq(dTr’~) = 0. 

Here, K q can be identified to a n x n block matrix, with blocks of size d x d, and the (i, j )—th 
bloc is given by e(a— Xj)Id , with Id the identity matrix. It is invertible (see[51]), and there exists 
S = a matrix of size n x n such that Y^k =1 S(xi,Xk)e(x k — Xj) = 1 if i = j 

and 0 if i ^ j. 

Then, for q = (xi,..., x n + m ), u = (iti,..., u n ) and p = (pi,... ,p n +m), we know that we have 

n+m 

K v CqP(x) = e(x-Xj)pf. 

3=1 

Hence, if (iq,..., v n ) = £, q K v £|p, we obtain 

n+m 

Vj = E e(xj -Xk)Pk, j = 1, • • •, n. 
k= 1 

So the equation K n ^u = ^ ^K v ^p becomes 

n n+m 

E e ( x 3 - x k)ul = E e ( x 3 - x k)pJ, i = 

k= 1 k= 1 

This linear equation is solved in u by 

n n+m n n+m 

Uj = EE S(xj,Xj)e(xj,x k )pk =Pi + E E S{xi,Xj)e{xj,x k )p k , 

j=l fc=l . 7=1 fc=n+l 

and therefore, if dpH^(q,p, u) = v £qU = (hi,..., v n +m), we have 


v a 


nun n+m 

^2e(x a ,Xi)uf =^2e(x a ,Xi)pl + ^2 E S{x i ,x j )e{x j ,x k )pl. 


2=1 2=1 2,.7=1 /c=n+l 

Noting that dqdn^ = 0, we must now compute 

dqH?(q,p,u) = (p- ^dir*-u){{dqKq)dir*-u) + (p - d7r|u) = (p - ^dn*-u)({dqKq)dTr*-u). 

If we denote p' = p - \dn*-u = ( p[ ,... ,p+ m ) and write -dqH*(q,p,u) T = (aq,... ,a n+m ), we 
have 




1 x ^ 

- E e(x Q - Xj)[Pa •« 1 +Uo- p'](x a - Xj) T 

a z ' 

2=1 

^ n+m 

- V e(x a — Xj)[u a • p'i](x a — Xi) T , a =1, ...,n, 

a . z ' 

2 =n+l 
1 n 

- E e(x a - Xi)ba • Ui](x a - Xi) T , a =n + 1,..., n + m. 

2=1 
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On the other hand, p' a = p a — i u a , hence, we obtain 

1 n 

a a = - e(x a - Xi)[p a ■ Ui + Ua ■ Pi - Ua ■ Ui\(x a - Xi ) 1 

rr ‘ ■* 


i=1 
n+m 


+ - V' e(x Q - Xi)[u 0 -pi](x Q - Xi) T , 

/T z - * 


i=n -\-1 


tt«= - V' e(x a - Xi)[p a • + (+ - Xi) T , 

(J Z -' 

i= 1 

We can finally write the lifted geodesic equations on S 

n n-\-m 


a =1,..., n, 


a =n + 1,..., n + m. 


Ua — Pa ^ ^ ^ $(x a: Xi)c(^Xi Xfc)pfc, 

i= 1 /c—n+1 
n 

y^e(x a - Xi)uf 
'■ 1 

n 

^ e(x Q — Xi) [p a ■ Ui + u a ■ Pi - u a ■ Ui] (x a - Xi) 1 


x a = 


a =1,..., n, 


a =1,..., n + m, 


Pa = 


2 = 1 

1 

(7 


1 


2=1 

n+m 


+ - V' e(x a - Xi)[u a -pi](x Q - x+ 

<T z ^ 


i=n+l 

l ” 

Pa ~ - e(x Q - Xi)[p a • lti](x a - Xj) T , 

(7 z ' 

2=1 


a =1,..., n, 

a =n + 1,..., n + m. 


Tangent spaces of landmarks with the Gaussian kernel In this example, we keep our 
Gaussian kernel K as above and S = Lmk n {S, d ) but we take S = TS = S x (R d ) n and 7r(g, v ) = q. 
We denote q = (q , v ) = (xi, ... ,x n ,vi, ... ,v n ). Again, we want to compute the geodesic equations, 
first by solving 

X ir(q)U +(<y)A" CqPi 

then by computing ddpH^q.p.u) = £qK v £*^u and dqH(q,p,u ) = (p — ^dTT*~u)((dqKq)dTT*~u). 
Again, the term (p — dTrtu)(Kq((dqdTT^)u) does not appear in dqH(q,p,u) because dq{dn*^) = 0. 

A momentum on S is denoted p = (p, l) = (pi,... ,p„, l \,..., l n ), and we have 

n 

Cq,v{pM x ) = y]pj(W(xi)) d+^dX+u,-)). 

f=i 

Therefore we get 

n i n 

KV £q,v(P, 00) = X+O - Xj)pj - - e(x - Xj)[(x - xj) ■ v+J. 

1=1 17 fe=i 

Hence, the equation K^^u = £|p becomes the linear problem 

n n i n 

^e(xj - x k )ul = ^e+- - x k )pj -^e(x^ - x fe )[(x i - x fe ) • p+J, j = l,...,n, 

k=l k= 1 ^ fc=l 
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solved by 


H =Pi -V' S(xi,xj)e(xj - x k )[(xj - x k ) ■ v k ]l k , i = l,...,n. 

nr • 


j,k= 1 


Then we compute t; q , v K l £*u = (w i,..., w n , w[,..., w ' n ) 

n 

Wa= y^ t e(x a - Xj)uf, 

2=1 

n i n 

= '£e(xa-x i )tf - ^2 e (Xj - Xk)[(Xj - x k ) ■ Vk]l 


2=1 


j,k=1 


w' a = -y e(x a - Xi)[v a ■ (x a - Xi)]uj, 

(7 z —' 


2=1 

Now if we denote 


a = 1,..., n. 


p' =P~ ^dn~U = = {Pi ~ ,...,p n ~ ^U n ,l U ...,l n ), 


we obtain 


p r (Kqd,Tr*u) = y e ( x i - Xj)\p[ ■ Uj] - "22 e i x i ~ x i)\ v i ’ ( x i - x o)][U ' u j\- 


i,j =i 


i,i=i 


Therefore, for -dqHT(q,p,u) = (ctq,..., a n , /?i,..., /?„), we get 


— 


y, e(x a - Xj)[Pa • Uj + Pj ■ U a }{x a - XjY 


0 =1 


1 


+ - y e(x Q - ^-)([Z a • Uj]Vg - • Uo\vJ) 


3 =1 


- 2 y e ( X a - 2r,-)[(x 0 - 2b) • ([ia • - [lj ' U a ]Vj)\{x a - Xj ) T , 


i=l 


/3 a — y ^ e(x a Xj) \lg ‘ 2/q] (x a Xj) 


a = 1,..., n. 


3 = 1 

Since p[ = Pi — \Ui , [p^ • + p' • u a ] = [p a • Uj + Pj ■ u a — Uj ■ u a \ and we obtain the geodesic 

equations 


1 ” 

ia=Pa - y S(x a ,Xj)e(Xj -Xk)[{Xj - Xk) -Vk]lk, 


j,k =1 


Xa = 


y e{x a - x i)uf , 


2=1 


y a =-y e(x a - Xj)[w a • (x a - Xi)]i 

nr • ■* 


a = 1,..., n, 
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and 


*| * It ' 

Pa = ~ yZe(x a - Xj)\p' a ■ Uj+p ' • u a ](x a - Xj) T 

3 = 1 

1 n 

+ - V e[x a - Xj)([la ■ Uj]vl - [lj ■ U a ]vj) 

a z ' 

3 =1 

x „ a = 1,..., n. 

~ ^2 - Xj) ■ {[la ■ Uj]v a - [lj ■ U a ]Vj)]{x a ~ Xjf , 

3 = 1 

1 n 

la — ^ ' c{Xa Xj ) [^a * Xj') , 

(T z ' 

i=i 

Conclusion 

This paper gave a comprehensive study of the geometry of abstract shape spaces in the LDDMM 
framework, and introduced several possible areas of applications, from fibered shapes to lifted 
shape spaces, to even more general shapes. We also showed that the LDDMM method is, at least 
on infinite dimensional shape space, largely a sub-Riemannian problem, not a Riemannian one. 

Other than those applications, several things still need to be investigated. First of all, the case 
of weak Riemannian (and even sub-Riemannian) structures, as in [13, 39]. Most of the results still 
work in this case, but they require additional attention, because the Hamiltonian equations do 
not always have a solution, so that the geodesic flow must be checked to exist on a case by case 
basis. The case of Frechet and convenient shape spaces was also only briefly mentionned, and may 
deserve additional attention. 

Another kind of problem left untouched is that of image matching. An image on M is a function 
I : M —> R. V s (M) acts by composition of the inverse on the right ip • I = I o <p~ l . This action 
does not satisfy the shape space hypothesis, and a different point of view must be adopted. Also 
see [16, 17] for another applications of sub-Riemannian geometry to image analysis. 

Finally, infinite dimensional sub-Riemannian manifolds are still very poorly understood, par¬ 
ticularly because of the appearance of elusive geodesics. 
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